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(U)  A theoretical  analyiit  of  additive  drag  of  two-dimensional  inleti 
hai  been  performed.  A pieliminary  review  of  tupertonic  inleta  provides  a 
description  of  the  various  types  of  inlets,  an  introduction  to  inlet  terminology 
and  performance  criteria,  and  the  definition  of  additive  drag.  The  application 
of  geometrical  streamline  tracing  is  introduced  for  simple  geometries  at  a 
method  of  determining  inlet  capture  area  ratio. 

0t)  The  spillage  from  the  sides  of  a two-dimensional  inlet  it  treated  by 
the  mothod  of  homogeneous  conical  flows.  This  method  Is  described  and 
extended  to  determine  the  three-dimensional  perturbation  velocities  within 
the  con(cal  fluwfleld. 

0s)  A computer  program  has  been  written  to  determine  the  shape  of 
the  streamtubo  captured  by  the  Inlot.  in  addition  to  tracing  the  streamlines, 
the  program  alto  numerically  integrates  the  pressures  over  the  streamtube  to 
calculate  the  additive  drag  with  tidetpillage. 

(0)  The  computer  program  hat  been  applied  to  aeveral  inlet  designs. 
Comparisons  of  spillage  are  made  with  experimental  data.  The  effects  of 
various  inlet  variables  are  considered.  These  variables  include  ramp  angle, 
Mach  number,  presence  of  sidewalls,  and  leading  edge  tweepback. 
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Chapter  1.  INTRODUCTION 

In  recent  years,  there  has  been  an  increasing  emphasis  in  supersonic  airbreathing 
propulsion.  Military  aircraft  are  being  designed  for  sustained  supersonic  flight.  There  is 
renewed  interest  in  ramjet  propulsion  for  missile  applications  at  supersonic  conditions. 
The  inlet  is  a major  component  of  these  airbreathing  propulsion  systems.  For  some 
applications,  two-dimensional  inlets  are  supplanting  the  standard  axisymmetric  inlet.  To 
determine  the  overall  performance  of  any  propulsion  system,  all  of  the  various  drag 
components  must  be  properly  accounted.  One  of  the  drag  components  attributed  to 
the  inlet  is  additive  drag.  Determination  of. additive  drag  is  relatively  simple  for  an 
axisymmetric  inlet  or  for  a simplified  model  of  a two-dimensional  inlet.  However, 
when  sidespillage  of  a two-dimensional  inlet  is  included,  the  determination  of  inlet 
mass  flow  and  additive  drag  becomes  complicated.  This  problem  is  considered  herein. 

Chapter  2 contains  a review  of  supersonic  inletr,  The  operation  and  terminology 
of  inlets  are  discussed,  and  additive  drag  Is  defined.  The  method  of  streamline  tracing 
to  determine  inlet  mass  flow  and  additive  drag  is  introduced.  Chapter  3 reviews 
previous  investigations  of  spillage  of  two-dimensional  inlets.  The  method  of 
homogeneous  conical  flow  is  introduced  and  extended  to  determine  the  perturbation 
velocity  components.  In  Chapter  4,  this  extended  theory  is  employed  in  a numerical 
technique  for  streamline  tracing.  A computer  program  is  described.  Applications  of  this 
computer  program  are  presented  in  Chapter  5.  Mass  flow  determined  by  this  method  is 
compared  to  measured  data.  The  program  is  applied  to  demonstrate  the  effects  of 
various  inlet  design  variables,  Most  of  the  conclusions  concerning  the  design  of  inlets  is 
included  in  this  chapter.  Chapter  6 summarizes  the  previous  work  and  dlscusscB  the 
range  of  applicability  of  the  computer  program. 


Chapter  2.  REVIEW  OF  SUPERSONIC  INLETS 


Before  proceeding  to  a discussion  of  spillage  in  two-dimensional  inlets,  it  would 
be  useful  to  review  the  operation  and  terminology  of  inlets  for  airbreathing  propulsion 
systems.  More  complete  discussions  are  presented  by  Faro  (Ref.  I)  and  the  staff  of 
Johns  Hopkins  Applied  Physics  Laboratory  (Ref.  2).  The  inlets  are  probably  the  most 
prominent  feature  distinguishing  an  airbreathing  propulsion  system  (such  as  a ramjet, 
turbojet,  or  alMugmented  rocket)  from  a rocket  propulsion  system.  Obviously,  the 
primary  function  of  an  inlet  is  to  capture  and  induct  air  as  a working  fluid  to  a 
desired  location  for  combustion  and  subsequent  expulsion  through  an  exhaust  nozzle. 
In  addition,  the  inlet  is  usually  required  to  decelerate  the  flow.  For  a turbojet,  the 
velocity  must  be  reduced  to  a value  compatible  with  the  compressor  blades,  a Mach 
number  of  about  0.5.  For  a ramjet,  the  velocity  of  the  air  must  be  reduced  to  the 
order  of  magnitude  of  the  flame  speed  or  to  the  order  of  a hundred  feet  per  second. 
In  a supersonic  inlet  this  is  achieved  by  a series  of  shock  waves. 
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Inlets  typically  are  of  two  types,  axiaymmetric  or  two-dimensional,  although 
special  circumstances  may  dictate  a more  complex  three-dimensional  inlet. 
Axisymmctnc  inlets  have  demonstrated  utility  and  applicability  in  some  of  the  earlier 
ramjet  programs  such  as  BOMARC  and  TALOS.  Axisymmetric  inlets  are  relatively 
inexpensive  to  build.  Structurally,  the  circular  duct  offers  the  best  shape  for  the 
pressure  loads.  This  permits  lighter  weight  construction.  The  boundary  layer  bleed 
requirements  for  an  axisymmetric  inlet  are  minimal.  However,  the  performance  of  an 
axisymmetric  inlet  decreases  with  angle-of-attack.  A two-dimensional  inlet  may  be 
oriented  to  achieve  improved  angle-of-attack  performance.  The  Mudto.  Ranrjet  Engine 
(MRE)  and  some  aircraft,  such  as  the  F14  and  FI5,  have  en  . .oyed  two-dimensional 
inlets  for  this  reason.  Two-dimensional  inlets  also  provide  packaging  flexibility  in  that 
any  combination  of  width  to  height  may  be  chosen  to  achieve  a given  inlet  area.  The 
two-dimensional  inlet  also  offers  versatility  in  the  choice  of  direction  of  compression. 
The  inlets  may  be  mounted  so  that  the  compression  turns  the  flow  inward  toward  the 
vehicle  or  outward  away  from  the  vehicle,  or  may  be  side-mounted  with 
downward-turning  Inlets.  The  compression  direction  may  be  selected  to  enhance 
angle-of-attack  performance  or  drag  characteristics.  For  variable  geometry  inlets,  the 
ramp  construction  permits  the  greatest  range  of  configurations  providing  near-optimum 
performance  over  u wide  flight  envelope.  The  two-dimonsionul  inlet  is  heavier  and  more 
expensive  to  fabricate  than  the  axisymmetric  inlet.  In  addition,  the  comer  flow  creates 
a thicker  boundary  layer  requiring  more  bleed. 

The  simplest  type  of  supersonic  inlet  is  a normal  shock  inlet  in  which  a normal 
shock  stands  at  the  front  of  the  Inlet,  as  shown  in  Figure  la.  With  subsonic  flow 
behind  the  shock,  the  (low  may  be  further  decelerated  by  a diverging  duct,  as  in  a 
subsonic  inlet.  This  ,nlet  is  typically  the  most  stable  and  least  sensitive  to  upstream 
conditions.  Satisfactory  performance  can  be  obtained  from  a normal  shock  inlet  at  low 
supersonic  speeds. 

The  diffuser  duct  may  be  converging-diverging  to  form  a reverse  de  Laval 
nozzle  us  is  shown  in  Figure  lb.  When  the  inlet  is  started,  the  supersonic  flow 
decelerates  in  the  converging  entrance  section,  passes  through  a normal  shock  just 
downstream  of  the  throat,  and  further  decelerates  subsonically  in  the  divergent  section. 
This  type  of  Inlet  has  a starting  problem.  In  the  unstorted  mode,  the  normal  shock 
will  stand  in  front  of  the  inlet.  The  throat  area  must  be  increased  or  the  inlet  must  be 
oversped  (If  possible)  in  order  to  "swallow”  the  shock. 

The  total  pressure  loss  of  a supersonic  inlet  n ay  be  reduced  by  means  of  an 
oblique  shock,  This  type  of  inlet,  shown  in  Figure  L may  be  either  a conical  inlet  or 
a two-dimensional  ramp  inlet;  many  of  the  conclusions  arc  applicable  to  both  kinds. 
As  indicated  in  the  figure,  an  initial  compression  surface  causes  an  oblique  shock 
which  decelerates  the  flow.  At  the  entrance  to  the  inlet  duct,  a normal  shock  occurs 
at  a Mach  number  lower  than  freeetream.  The  product  of  the  total  pressure  recovery 
of  the  oblique  and  normal  shocks  is  greater  than  that  of  a single  normal  shock  at 
freestream  conditions.  Figure  2 indicates  the  total  pressure  recovery  of  conical  Inteti 
for  several  Mach  numbers.  As  can  be  seen,  there  is  an  optimum  cone  angle  for  each 
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Mach  number.  This  cone  angle  may  be  presented  as  a function  of  Mach  number  as  in 
Figure  3.  The  pressure  recovery  of  the  oblique  shock  inlet  may  also  be  increased  by 
internal  contraction  which  reduces  the  Mach  number  prior  to  the  normal  shock  at  the 
throat  section.  Internal  contraction  it  limited  by  inlet  starting  considerations.  Finally,  a 
subsonic  diffuser  section  provides  further  deceleration  of  the  flow, 

Since  the  single  oblique  shock  improved  the  performance  of  the  inlet,  it  may 
be  suspected  that  an  additional  oblique  shock  before  the  normal  shock  may  also  be 
beneficial.  This  is  indeed  the  case.  A double  cone  or  ramp  inlet  is  illustrated  in  Figure 
Id.  As  in  the  single  oblique  shock  inlet,  there  is  an  optimum  combination  of 
deflection  angles.  Figure  4 presents  the  total  pressure  recoveries  of  various 
combinations  of  cone  angles  for  a double  cone  inlet.  The  optimum  cone  angles  are 
shown  as  a function  of  Mach  number  In  Figure  5, 

Continued  improvement  in  total  pressure  recovery  may  be  obtained  by 
increasing  the  number  of  oblique  shocks.  Figure  6 demonstrates  the  total  pressure 
recovery  for  optimum  combinations  of  cone  angles  for  multiple  oblique  shocks.  The 
number  of  shocks,  n,  includes  the  terminal  normal  shock  in  addition  to  the  oblique 
shocks.  In  the  limit,  the  oblique  shocks  become  an  isentroptc  compression  shown  In 
Figure  Ic.  Since  a true  isentropic  compression  requires  an  initial  zero  slope,  the 
Isentropic  inlet  starts  with  an  initial  nomzero  slope,  either  a selected  angle  or  the  angle 
which  yields  a selected  pressure  recovery  at  design  conditions.  The  coordinates  of  such 
Isentropic  surfaces  have  been  calculated  for  various  design  Mach  numbers,  An  example 
for  conicul  Inlets  is  included  as  Figure  7.  In  this  figure,  focal  points  are  indicated  for 
each  design  Mach  number.  At  the  design  Mach  number,  the  conical  shock  from  the  tip 
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and  all  compression  waves  generated  front  the  surface  coalesce  at  this  point.  In 
practice,  this  is  usually  made  to  be  the  location  of  the  cowl  lip.  All  dimensions  in  the 
figure  are  ratioed  to  the  radius  of  the  focal  point  Yp  so  that  when  the  cowl  diameter 
is  selected,  the  inlet  shape  is  specified. 

The  relative  performance  of  air  inlets  is  compared  on  the  basis  of  three  primary 
figures  of  merit.  These  arc  (I)  capture  area  ratio,  (2)  total  pressure  recovery,  and  (3) 
drag. 


The  capture  area  ratio  indicates  how,  much  air  is  delivered  to  the  propulsion 
system.  The  mass  flow  at  any  point  may  he  written  as  pAV  where  p is  the  density,  A 
is  the  cross-sectional  area,  and  V is  the  velocity.  When  referenced  to  frccstream  or  0 
(zero)  conditions  where  the  density,  p0,  and  velocity,  V0,  are  readily  obtained,  the 
mass  flow  (which  may  be  measured  elsewhere)  defines  the  freestream  area  of  the 
captured  streamtube,  A«.  The  captured  muss  flow  may  be  normalized  by  dividing  by 
the  mass  flow  that  could  be  delivered  by  the  area  of  the  inlet  or  cowl  area,  Ac.  This 
mass  flow  is  P0ACV0.  Thus,  the  ratio  of  the  uctual  mass  flow  to  the  maximum 
possible  becomes  the  capture  area  ratio,  A0/Ac.  In  addition  to  the  physical  meaning, 
based  upon  a measured  mass  flow,  the  capture  area  may  be  determined  theoretically 
by  geometrical  analysis.  Consider  a single  ramp,  two-dimensionul  inlet  as  shown  In 
Figure  8.  The  cowl  area,  A , is  simply  the  cowl  height,  h , multiplied  by  the  inlet 
width,  W.  The  flow  behind  the  shock  is  directed  parallel  to  the  compression  surface  at 
an  angle,  0%,  A streamline  may  be  drawn  from  the  cowl  lip,  parallel  to  the 
compression  surface,  to  the  intersection  with  the  oblique  shock  inclined  at  wave  nngle, 
0W.  Upstream  of  the  shock  wave,  the  flow  is  horizontal.  Thus,  the  captured 
streamtube  of  height,  h0,  and  width,  W,  for  two-dimensional  Row  is  defined.  The 
cupture  area  ratio,  A0/Atf,  is  equal  to  the  ratio  of  heights,  U0/he,  in  this  case. 
Extension  to  multiple  rainps  requires  only  identification  of  the  flow  direction  behind 
each  shock  and  tracing  the  streamline  from  the  cowl  back  to  freestream.  An 
axisy mmetric  inlet  is  only  slightly  more  complicated  in  that  the  flow  is  not  parnllel  to 
the  surfuce  but  vuries  between  rays.  The  streamline  must  be  traced  from  ray  to  ray  to 
the  initiul  shock.  Following  this  principle,  the  capture  area  ratio  of  any  inlet,  even  a 
three-dimensional  inlet,  may  be  determined  by  a streamtube  truce  if  the  flowfleld  is 
known.  Mass  removal,  such  as  by  sidespill  or  boundary  layer  bleed,  may  also  be 
referenced  to  the  possible  mass  flow  to  yield  an  area  ratio  for  each  mass  removal 
process. 

The  second  figure  of  merit  has  been  alluded  to  in  the  description  of  the  inlets. 
The  total  pressure  recovery  is  a direct  measure  of  the  maximum  pressure  which  can  be 
achieved  in  the  combustion  chamber  and  hence  of  the  thrust  developed. 

There  are  nine  possible  modes  of  operation  of  an  inlet  resulting  from 
combinations  of  three  possibilities  for  each  of  two  independent  variables.  These  nine 
modes  are  illustrated  in  Figure  9.  First,  an  inlet  may  operate  above,  at,  or  below  the 
design  Mach  number.  This  is  presented  in  the  vertical  sequence  in  the  figure.  Typically, 
the  design  Mach  number  Is  the  Mach  number  at  which  the  oblique  shock  of  the  inlet 
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intersects  the  cowl  lip.  For  multiple  shock  inlets,  the  ramp  contour  is  designed  so  that 
all  shocks  coalesce  at  the  cowl  Up.  Since  the  shock  angle  for  constant  deflection 
decreases  with  increasing  Mach  number,  the  shock  passes  above  the  cowl  below  the 
design  Mach  number  and  enters  the  cowl  above  design  Mach  number.  The  effect  of 
Mach  number  on  total  pressure  recovery  and  capture  area  ratio  are  illustrated  on  the 
right  side  of  the  figure.  Total  pressure  recovery  typically  decreases  with  increasing 
Mach  number,  slowly  below  design  and  more  rapidly  above.  Often  a sharp  break  in  the 
curve  may  be  noted  at  the  design  Mach  number.  Below  design,  analysis  of  a 
two-dimensional  inlet  is  straightforward  in  the  manner  previously  described.  Above 
design,  analysis  is  complicated  by  the  non-uniform  conditions  at  the  cowl  lip  plane. 
Part  of  the  flow  has  passed  through  the  oblique  shock  system  with  a loss  of  total 
pressure  and  velocity  and  a change  in  direction.  The  remainder  of  the  flow  at  the  cowl 
Up  plane  is  still  in  the  undeflected  freestream  condition.  An  analytical  model  should 
include  some  averaging  technique  to  obtain  mean  conditions  at  the  cowl  lip  plane  for 
gross  inlet  behavior.  The  capture  area  increases  with  Mach  number  below  design.  Above 
design  Mach  number,  conditions  at  the  cowl  lip  ore  freestream  conditions.  Since  the 
flow  is  not  deflected,  the  area  of  the  captured  streamtube  is  the  same  as  the  inlet  area 
and  the  capture  area  theoretically  equals  1.0.  The  change  in  operation  appears  as  a 
sharp  break  in  the  capture  area  curve  at  the  design  Mach  number. 

The  other  variable  of  inlet  operation  is  whether  the  inlet  is  critical, 

supercritical,  or  subcritical.  These  conditions  are  illustrated  in  the  horizontal  sequence 
of  Figure  9.  In  the  critical  condition,  the  terminal  nonnal  shock  is  at  the  cowl  lip 

plane  of  an  external  compression  inlet  or  at  the  throat  of  a mixed  external-internal 
compression  inlet.  In  the  supercritical  condition,  the  terminal  normal  shock  is 

downstream  of  the  critical  position,  occurring  in  the  subsonic  diffuser.  Due  to 

supersonic  expansion,  the  terminal  shock  occurs  at  a higher  Mach  number  than  at  the 
throat,  and  the  total  pressure  recovery  is  less.  Sinco  the  external  shock  system  is 
unchanged,  the  capture  area  ratio  remains  the  same  as  at  critical.  In  the  subcritical 
condition,  the  terminal  normal  shock  occurs  upstream  of  the  critical  position,  on  the 
compression  surface.  This  shock  location  results  in  a reduction  in  mass  flow  or  capture 
area  ratio.  Depending  upon  the  design  of  the  compression  surface,  subcritical  total 
pressure  recovery  may  be  greater  than,  equal  to,  or  less  than  critical  recovery.  The 
sequence  is  summarized  in  the  graph  at  the  bottom  of  Figure  9,  with  I representing  a 
subcritical  condition,  2 critical,  and  3 supercritical.  As  may  be  observed  from  the 
figure,  the  critical  condition  may  be  defined  as  the  condition  of  maximum  pressure 
recovery  at  the  maximum  capture  area  rutio.  The  location  of  the  terminal  nonnal 
shock  is  determined  by  the  combustor  temperature  and  therefore  fuel  flow  to  the 
combustor.  With  increased  temperature,  a higher  total  pressure  is  required  for  the  same 
mass  flow  rate.  To  achieve  higher  total  pressure,  the  shock  moves  toward  the  critical 
position.  For  further  combustor  temperature  Increases,  the  inlet  delivers  less  mats  by 
subcritical  operation,  thus  reducing  the  pressure  recovery  requirements. 

The  primary  concern  in  testing  inlets  is  the  determination  of  the  critical 
condition.  It  is  generally  assumed  that  the  inlet  may  be  easily  operated  supercritically 
for  less  total  pressure  recovery.  In  subcritical  operation,  an  instability  known  as 
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“diffuser  buzz"  may  occur.  During  buzz  the  shock  position  is  unstable,  oscillating  from 
supercritical  to  subcritical  operation,  possibly  to  a detached  position.  This  oscillation  is 
accompanied  by  static  pressure  oscillations  throughout  the  propulsion  system.  To  avoid 
the  onset  of  buzz,  propulsion  systems  are  designed  to  avoid  the  unstable  region.  For 
aircraft  turbojet  engines,  where  a variable  geometry  inlet  is  economically  feasible,  the 
inlet  configuration  is  adjusted  at  each  flight  condition  to  maintain  near  critical  or 
slightly  subcritical  operation  to  obtain  maximum  pressure  recovery.  This  may  be 
achieved  by  change  of  ramp  angles  for  a two-dimensional  inlet,  by  change  of 
centerbody  location  for  an  axisymmetric  Inlet,  or  by  varying  the  bleed  flow. 
Subcritical  operation  has  the  additional  advantage  of  providing  relatively  undistorted 
flow,  desirable  for  a turbojet  compresaor.  For  a fixed-geometry  inlet,  typical  of  ramjet 
application,  the  Inlet  is  operated  supercritlcally.  The  inlet  is  designed  with  a pressure 
margin  or  maximum  pressure  recovery  less  than  critical.  This  acts  as  a safety  factor  to 
assure  that  the  inlet  will  not  achieve  critical  condition,  much  less  an  unstable 
subcritical  condition. 

The  third  figure  of  merit  la  the  performance  cost  of  the  Inlet,  the  drag.  There 
are  several  types  of  drag  attributed  to  the  inlet.  There  is  wave  drag  associated  with  the 
cowl  shape  and  on  the  cowl  lip.  For  two-dimensional  inlets  with  sideplates,  there  is  a 
corresponding  sidcplate  lip  and  wave  drag.  For  inlets  Installed  on  the  aft  portion  of  a 
vehicle,  a boundary  layer  diverter  It  provided  to  raise  the  inlet  above  the  vehicle 
boundary  layer  so  that  it  does  not  ingest  the  low  pressure  recovery  air.  This  is  usually 
a wedge-shaped  structure.  There  Is  a drag  associated  with  the  diverter;  and  an  optimum 
diverter  height  may  be  found  by  trading  off  pressure  recovery  gains  with  Increased 
drag.  Mass  removal  processes  of  bleed  and  sidespill  also  have  momentum  losses  and 
hence  drag  associated  with  them.  Finally,  there  is  the  additive  drag  of  an  inlet. 

The  term  additive  drag  was  originally  defined  by  Ferri  and  Nucci  (Ref.  3).  A 
mathematical  description  of  additive  drag  hat  been  presented  by  Sibulkin  (Ref.  4). 
First,  consider  a simple  ramjet  engine  in  accelerated  flight  at  zero  angle-of-attack,  as 
shown  in  Figure  10a.  The  net  propulsive  thrust  is  the  resultant  sum  of  the  axial 
components  of  the  pressure  and  friction  forces.  The  net  internal  thrust,  Fn  ,,  may  be 
evaluated  from  the  momentum  change  between  the  entrance,  Station  1,  and  the  exit, 
Station  e. 


n.i 


v*. 


(2-1) 


where  F j ■ stream  thrust  at  j « mVj  + Aj(p,  - p0).  Thus,  the  net  propulsive  thrust, 
Fp,  Is  equal  to  the  net  internal  thrust  minus  tne  external  drag  forces,  Fd. 


n.i 


(2-2) 


In  practice,  engine  performance  Is  usually  determined  from  freestream  and  exit 
conditions  instead  of  from  inlet  entrance  conditions.  Thus,  referring  to  Figure  10b, 
which  includes  the  streamtube  from  0 to  1,  the  net  thrust,  Fn,  due  to  momentum 
change  between  Stations  0 and  c is: 
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In  this  case,  however, 


(2-3) 


since  the  momentum  change  between  Stations  0 and  1 has  not  been  considered.  This 
momentum  change  may  be  considered  at  an  additional  drag  force,  D|(  called  additive 
drag.  Thus, 


c a p _ p _ n 
rp  rn  rd  • 

which  defines  additive  drag.  Comparing  Equations  2-1  and  24  yields 

D.-Fn-Fn(t-/Wo 


(24) 


(2-5) 


or 


D,-mV,  +A,(p,  - p0)-  mVQ 

Considering  the  streamtube  from  0 to  1 as  indicated  in  Figure  10c, 
/•Sta  1 

/St  0 (P  " p0)dA*  “ Al(pl  “ P0)"-mV0  + mV, 

where  dAx  ■ axial  projection  of  surface  area,  or  from  Equation  2-6, 
•Sta  1 


/•sta  l 

D.“  / (P“P0>dAs 
J Sta  0 


(2-6) 


(2-7) 


(2-8) 


This  integral  form  of  the  equation  is  applicable  to  other  geometries.  Evaluation 
of  the  momentum  terms  results  in  a slightly  different  form  of  additive  drag  for  each 
situation.  Mount  (Kef.  5)  presents  the  following  equation  for  both  the  axisymmetric 
and  two-dimensional  oblique  shock  inlets  shown  i.  Figure  1 1 . 

/■Sta  1 

D«"7c*  n (P-Po*^*- 

/ StaO 

- gV0  + (V  Po>A.  («> 

where  p(  is  the  average  pressure  acting  on  the  compression  surface,  A,. 


w 


gV,  +(p,  - Pq)A|  j cos  9 


.0 


f v,  cot 9 


bt-P0»A,  cote 


~ V,  coif 
to,  -PQ>  A,  coif 


■#,-P0>\ 
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This  differs  from  the  relation  for  the  open  nose  inlet,  Equation  2-6,  by  the 
addition  of  the  term  (fi,  - p0)A(1  which  is  the  pressure  force  on  the  compression 
surface,  and  by  a correction  term  of  cos  6 to  the  exit  momentum  to  account  for  a 
m .' -axial  exit, 

Additive  drag  appears  to  be  only  a bookkeeping  term.  Its  physical  significance 
may  be  seen  by  considering  the  boundary  of  the  captured  streamtobe  to  be  replaced 


J 'U t+L* >-•  - > n..-, 
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by  a functionless  sheet  of  solid  material.  No  flow  may  pass  through  the  sheet  and 
since  it  conforms  to  the  shape  of  the  streamtube,  the  flowfleld  remains  undistubed. 
The  pressure  forces  acting  on  this  extended  cowl  are  easily  visualized.  Removal  of  this 
extension  does  not  change  the  flow  and,  hence,  the  momentum  through  the  inlet 
system.  Distortion  of  the  cowl  extension  would  alter  the  flowfleld  and  result  in  a 
different  pressure  distribution. 

Additive  drag  occurs  in  two  situations.  When  the  inlet  is  operating  critically  or 
supercritically,  at  or  above  the  design  Mach  number,  the  leading  edge  shock  is  at  or 
within  the  cowl  lip.  Up  to  the  cowl  Up,  the  .How  is  undisturbed.  Hence,  the  flow  area 
does  not  change,  and  the  additive  drag  is  zero  for  purely  two-dimensional  flow.  With 
suberitieu!  operation,  the  terminal  normal  shock  occurs  on  the  compression  surface  and 
subsonic  flow  occurs  behind  the  shock,  Osmon  (Ref.  6)  has  considered  subcritical 
additive  drag  In  the  transonic  region,  He  employed  Moeckel’s  (Ref.  7)  method  for 
prediction  of  the  plane  shock  shape.  This  equation  was  modified  to  account  for  lateral 
spillage;  various  coefficients  were  adjusted  to  fit  experimental  data,  With  a 
determination  of  the  terminal  shock  shape  and  position,  the  pressures  behind  the 
various  shocks  were  used  to  determine  subcritical  udditivc  drag. 

The  second  condition  of  additive  drag  is  supercritical  additive  drag.  In  this  case, 
the  inlet  Is  operating  critically  or  supercritically  below  the  design  Much  number. 
Supercritical  additive  drug  can  be  evaluated  for  simple  geometries  by  meuns  of  the 
pressure  integral.  Consider  the  single  ramp  two-dimensional  Inlet  of  Figure  8.  The  static 
pressure  behind  the  shock  Is  constant  along  the  streamline  and  may  be  obtulned  from 
plane  oblique  shock  relationships.  The  area  over  which  this  pressure  acts  is  simply 
(Ac-  Ap)  or  Acd  - A0/Ai;).  For  multiple  oblique  shocks,  this  process  may  be  repeated 
for  each  shock.  For  axisymmetric  Inlets,  the  same  approueh  is  applied,  numerically 
integrating  acioss  the  rays.  Sharp  and  Howe  (Ref.  8)  have  employed  this  technique  for 
determination  of  additive  drag  In  an  automated  procedure  for  analysis  of  inlets.  For 
two-dimensional  inlets  of  finite  width,  there  is  some  lateral  spillage.  Sideplatcs  may  be 
employed  to  restrict  or  climinntc  this  lateral  flow  at  the  expense  of  Increased  weight 
and  drug.  Sidewalls  may  also  be  constructed  to  provide  additional  sidewise 
compression,  If  a sidcplate  is  designed  to  contain  the  entire  shock  system,  there  will 
be  no  sldespill.  Sideplatcs  urc  typically  constructed  with  the  sidcplate  leading  edge  ns  a 
straight  line  from  the  cowl  Up  to  the  inlet  leading  edge.  Above  design  Mach  number, 
there  will  be  no  spillage  from  this  type  of  sidcplate.  In  general,  however,  there  will  be 
lateral  spillage  from  a two-dimensional  inlet.  This  will  alter  the  flowfleld  and  the  shape 
of  the  captured  streamtube.  To  determine  additive  drug,  the  shape  of  the  enptured 
streamtube  must  be  determined.  This  requires  knowledge  of  the  flowfleld  about  the 
compression  surface  with  lateral  spillage. 
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Chapter  3,  PREVIOUS  INVESTIGATIONS 

Little  theoretical  work  has  been  performed  on  the  problem  of  sidespill  in 
two-dimensional  inlets.  This  is  partially  due  to  the  relatively  recent  emphasis  on 
two-dimensional  inlets;  most  of  the  previous  inlet  studies  had  been  directed  toward 
axisymmetric  inlets,  In  addition,  most  inlet  investigations  are  primarily  experimental. 
Preliminary  performance  estimates  are  determined  by  a combination  of  theoretical  and 
empircal  techniques.  Wind  tunnel  tests  are  then  performed  to  verify  the  performance 
estimates  and  to  select  various  design  features  such  as  bleed  system,  design,  boundary 
layer  diverter  height,  and  installation  effects. 


INVESTIGATION  OF  CARRIERE  AND  LEYNAERT 

Curriere  und  Leynuert  (Ref,  9)  have  considered  the  problem  of  mass  flow  in 
two-dimensional  inlets.  No  attempt  has  been  made  to  determine  the  effects  on  additive 
drag.  Consider  a single  ramp  inlet  with  ruinp  angle,  'if,  as  shown  in  Figure  12.  The 


FIGURE  12.  TWDtansnstonal  Inlet  Nomeudatuia. 
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inlet  hag  two  lateral  sidewalls  with  leading  edges,  AB  and  A'B\  at  an  angle  $.  It  is 
assumed  that  the  ramp  angle,  is  small  and  that  the  shock  waves  are  comparable  to 
Mach  waves.  This  permits  the  use  of  linearized  theory. 

Three  possible  situations  may  be  ‘ ssumed  a priori  depending  upon  the  value  of 
the  parameter  0W/£,  where  0 is  the  cotangent  of  the  Mach  angle,  m,  end  ia  equal  to 

n/m*-  1.  The  parameter  0W/S  characterizes  the  lateral  aspect  ratio  of  the  inlet  and  is 
the  ratio  of  the  width  of  the  inlet  to  the  radius  of  the  Mach  cone  at  the  inlet 
entrance  station.  The  three  possibilities  arc  illustrated  in  Figure  13.  The  Mach  cones 
may  be  completely  separate  or  may  overlap.  If  the  overlap  occurs  within  the 
boundaries  of  the  inlet  or  1 < 0W/8  < 2,  the  sldespill  is  unaffected  and  remains  the 
same  as  the  separated  case.  For  narrow  inlets  where  0W/8<  1,  the  Mach  cone 
generated  by  one  edge  influences  the  flow  at  the  opposite  edge, 

Consider  a uniform  approach  flow  with  velocity  components: 

u,-U0 
u,  ■ o 

u3«0  (3-1) 

After  a small  perturbation  such  as  that  imposed  by  the  inlet  ramp,  the  velocity 
components  are 


U|  * U0  + u' 


u2  « v 
u3«w' 


(3-2) 


or  normalizing  by  U0' 


u,  ■ U0(l  + u) 
u2  - U0v 
u3  •*  UQw 


(3-3) 


In  small  perturbation  theory,  it  is  assumed  that  the  perturbation  component! 
are  small  or 

u,  v,  w « 1 

Considering  the  inlet  in  Figure  12,  the  flow  enten  the  inlet  through  the  region 
BCB'C'  of  surface  A - W(h  - i/ffi)  by  geometry.  The  maaa  flow,  m,  after  linearization  by 
small  perturbations  is 


m ■ JJf>(  1 + u)U0dA 


(3-4) 


where  p is  the  mass  density  and  dA  an  element  of  the  control  surface. 
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The  local  density,  p,  may  be  related  to  freestream  density,  p0>  by  Mach 
number  functions: 


(3.5) 

p0  (1  + *-§-*■  M2)1^"1 

The  local  Mach  number,  M,  may  be  written  in  terms  of  perturbation  velocities 
M2  * |(U0  + u')2  + v'2  + w'2  ] /a2  (3-6) 

M2  ■ (Ug  + 2U0u')/a2  to  first  order 
the  speeds  of  sound  are  related  by  temperature  ratios 


a2  1RT  _ T _ 1 +1i  * Mo 

Substituting  a2  into  Equation  3-6 

Uj  +2U0u'  l +^M2 


M2 


“0  l+^Mj 


1 +3Lj^M2 


M2  ■ Mq(1  + 2u)  ZTT'  "T 
1 


Solving  for  M2 


M2(l+W 
M ■ 


Thus, 


1 - 2urjpMj 


y — i , * +rr-Mo 

+ W- — *rr, 

* i - 2u24-lmJ 


(3-7) 


(3-3) 


(3-9) 
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or 


1 +rr'Mol  „ 7-  iw2 

— - 2UV«0 

From  Equation  3-5 

£-<l-2u^M  . 

Expanding  with  a binomial  expansion  to  first  order: 

■ 1 - uM„ 

Pq  0 

Thus  Equation  3-4  may  be  written 

m -pQU0  //'  + u)(l  -uMg)dA 


(3-10) 


(3-11) 


(3-12) 


(3-13) 


m»pQU0  J^[  \ +u(l  - M§)]dA  to  first  order 
Since  -1,  Equation  3-13  becomes 


(3-14) 


Thus,  the  flow  rate  entering  the  inlet  may  be  determined  once  the  distribution 
of  the  axial  perturbation  velocity,  u,  is  found. 

Consider  the  pure  two-dimensional  or  00  case.  For  a weak  two-dimensional  wave 
where  the  flow  is  deflected  through  angle 


(3-15) 


* Refer  to  Section  4-7,  "Week  Oblique  Shockn,"  (n  Ehnumt  of  Gaidynamlct,  by  H.  W.  Llepmum  end  A. 
Rothko.  New  York,  N.Y.,  John  WUoy  A Som,  Inc.,  1962,  pp.  92-93. 
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Substituting  into  liquation  3-14 


I 


t; 

s 


i 


m~«p0U0  ffl  » + Wo>dA 

A 

m"  «p0U0(l  + 00O)A  (3-16) 

The  entrance  area,  A,  has  been  determined  by  geometry  to  be 

A = W(h-ipt)  (3-17) 

Since  h/C  =*  tan  { 

A - Wh(l  - \J/  cot  {)  (3-18) 

Thus  Equution  3-16  becomes 

tn°°  - p0U0Wh[  I - *(cot  t-  /}„)-  *2p0  cot  t)  (3-19) 

Carrie  re  and  Leynuert  present  this  equation  only  to  first  order: 

m"  - p0li0Wh[  1 - »Kcot  * - (30))  (3-20) 


The  ignored  term,  however,  is  not  negligible.  Consider,  for  example,  the  case  of 
a 10-degree  ramp  and  a 19.6-degree  sidewall  at  a freostream  Mach  number  of  1.96. 
The  term  4i2IS0  cot  f is  equal  to  0.144  which  is  not  negligible  as  compared  to  the 
other  terms.  T.he  value  of  nY“/p0U0Wh  is  determined  to  be  0.660  using  Equation  3-19 
and  0.804  using  Equation  3-20.  Thus,  neglecting  the  second  order  term  appears  to  be 
incorrect. 


Carrien:  and  Leynuert  present  ull  results  relative  to  the  two-dimensional  result, 
m"°.  Calculations  by  Fenain  by  the  method  of  homogeneous  flows  are  cited  without 
reference;  and  the  following  results  are  presented  without  proof. 


a.  With  lateral  plates  and  0oW/K  > 1 


m-  m 

1r~ 


\ l>  R k2K'2  - E 


>2 


\^r 


(3-21) 
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where 

■i  = P0  tan  J 

k'-v/T^P 

K',  H' d Complete  elliptic  integrals  of  the  first  and  second  kind,  respectively 
of  modulus  k# 

b.  Without  lateral  plates  and: 

1:  ff0w/g » k,  > i 


2.  k*  < k)  < 1 


m - m"* 
rii- 


ifi 

W1 

. 


n-k,)2» 

— t — 7 


3. 


k,  < k’  < I 


m-  rn 


' *W 


k’+as^. 


/7r— ? <1+kl>2 
- y/k  - k? £ sin 


-I 


*7 


, ki  2ki  - 
+ 2kj  cos" 1 jjr  + ~y~ 


(3-22) 


(3-23) 


(3-24) 


Note  that  for  a selected  inlet  height  and  Mach  number,  all  of  these  results  have 
the  form  m - tvF/wT  ■ constant  times  (\ph rM8/W)  which  is  expected  for  linearized 
theory. 


Application  of  these  equations  has  revealed  that  the  equations  cannot  be 
evaluated  when  the  cowl  angle,  £,  is  greater  than  the  Mach  angle,  p.  Physically,  this 
would  occur  when  the  Mach  wave  from  the  leading  edge  enters  the  cowl.  This 
observation  suggests  that  no  distinction  has  been  made  between  a Mach  wave  angle  and 
the  actual  shock  wave.  This  is  considered  to  be  in  error  since  there  is  a significant 
difference  between  the  Much  wave  angle  and  the  angle  of  a shock  wave  even  for  small 
ramp  angles. 
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HOMOGENEOUS  CONICAL  FLOW 

The  method  of  homogeneous  conical  flows  may  be  applied  to  the  problem  of 
spillage  of  two-dimensional  inlets.  Pawlikowski  (Ref.  10)  has  employed  this  method  to 
determine  perturbation  velocities  and  spillage  on  the  compression  ramp  of  the  inlet. 
Milne-Thomson  (Ref.  11)  presents  a basic  discussion  of  conical  flow.  He  employs  this 
method  to  determine  perturbation  velocities  on  various  airfoils.  Since  Milne-Thomson  is 
primarily  concerned  with  aerodynamic  forces,  he  uses  these  velocities  to  determine 
lifting  pressures.  The  analysis  of  Milne-Thomson  will  be  described  here  as  an  aid  to 
understanding  and  as  a contrast  to.  the  investigation  of  Pawlikowski.  Later,  the  velocity 
solution  of  Milne-Thomson  will  be  applied  to  the  spillage  problem. 

Conical  flow  is  a flow  pattern  in  which  the  state  of  the  flow  is  the  same  at 
every  point  along  a ray  drawn  from  an  origin  or  vertex.  Thus,  conditions  depend  only 
on  the  direction  of  the  ray,  but  not  on  the  distance  along  it.  Tills  means  that  any 
property  is  a homogeneous  function  of  degree  zero  in  x,  y,  z,  and  depends  only  on 
the  ratio  of  x:y:z.  A conical  boundary  condition  implies  conical  flow.  Since  the  edge 
of  an  inlet  compression  ramp  imposes  a conical  boundary  condition,  the  conical  flow 
analysis  may  be  applied  to  the  inlet  flowfield. 

For  small  perturbations  in  the  coordinates  of  Figure  14,  u,  v,  w < V,  and 
neglecting  squares  and  products,  the  perturbation  equation  for  supersonic  flow  is 


-(M2,  - I) 


•is . . d^s-rt 

3x2  3y2  3z2  ° 


(3-25) 


where  S may  be  any  of  the  quantities:  velocity  potential,  4>,  perturbation  velocity 
components  u,  v,  or  w,  pressure,  p,  or  density,  p. 


Equation  3-25  may  be  rewritten  as 


2 32S  d2S  32S  „ 

- cot2  p — T + — f + — x ■ 0 
Ox2  3y2  3z2 

since  cot  ju  " >/m2  - 1 

Substituting  X ■ ix  tan  a>,  Equation  3-26  reduces  to  the  Laplace  equation 
32S  d2S  32S  . 

— r + — r + — r « 0 

3X2  dy2  3z2 


In  spherical  polar  coordinates  as  defined  in  Figure  15, 

X ■ Rj  cos  w 
y * Rj  cos  0 sin  cj 
z * R(  sin  0 sin  u> 


(3-26) 


(3-27) 
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Laplace's  equation  becomes 


3 3S\.  1 3 / , 3 1 d2S_„ 

5K7(r.  5iqj+airs5«r 


(3-28) 


When  S is  independent  of  R, , aa  in  conical  flow,  Equation  3*28  becomes 


By  substituting 


Xs  log  tan  l/2w 

as  as  _ i ds 

3w  ■ cosec  w 5J  ifiTw  3x 
Equation  3*29  becomes 

ai|  + ^-o 

0X2  3fl2 

This  has  the  general  real-valued  solution* 

S * F(x  + 10)  + F(x~  i«) 

where  F is  an  arbitrary  twice  differentiable  function. 

ttn  1/2w,r+cww 


eX+i»  . 0xeiO  - c‘9  tan  l/2w  - (cos  0 + i sin  d)('fySwS>) 


X+W  m cos  0 sin  cj  + i sin  0 sin  u 
1 " I + cos  CO 


(3-29) 


(3-30) 


(3-31) 


(3-32) 


(3-33) 


from  definition  of  polar  coordinates 


y + iz 

trr; 


(3-34) 


* Th»  moil  $e>uenl  toluUon  li  S ■ F(x  ♦ i#)  + 0(x  - **)  but  Um  moit  asMiti  rttl  Mhitd  tolutioA  u 

Slvin  in  Equation  3-32, 


'•*‘*’*»*t*"WW*^*1 
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Thus  x + i0  is  a function  of  (y  + iz)/(X  + R,)  and  from  Equation  3-32,  the 
general  real-valued  solution  of  Laplace’s  equation  of  zero  degree  is  SQ  where 


(3-35) 


where  fj  denotes  an  arbitrary  function.  Since  R^  * X2  + y2  + i?  or 
R,  •s/y2  +z2  - x2  tan2  p ' 

X + R>  “ ix  tan  n +v/y2  + z2  - x2  tan2p 


where 


X + R,  ■ i(x  tan  n + R2) 


■ x2  tan2/i-  y2  - z2 


(3-36) 


X±iz  .f_n_ JLtiL. 
xVR,  1 ux  tan  u f 


u + r2 


(3-37) 


The  genera)  real-valued  solution  of  degree  zero  of  Equation  3-25  Is 


* \x  tan  n + Rj/  \x  tan  \i  + 

whero  f is  an  arbitrary  twice  differentiable  function.  This  solution  applies  to  all  points 
such  that  Kj  is  positive.  This  is  all  points  within  the  Mach  cone 


y2  + z2  ■ x2  tan2p 


(3-39) 


Introducing  Busemann'i  (Ref,  12)  complex  variable,  e*xtanjU  + R2,  the 
general  solution  becomes 

S0-f(«)  + T(3)  (340) 
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Introducing  the  dimensionless  variables,  r,  s,  s*,  which  satisfy 


then 


where 


s ■ r cos  0 ■ jnjjjTjJ 


s + ls*  rtcos  B + 1 sin  fl)  n m 

- ,2 -,*2)1/2  | +<1-12)1/2 


(3-41) 


(342) 


(343) 


If  the  point  re19  is  marked  in  an  Argond  diagram  called  the  physical  plane  and 
the  point  Re19  marked  in  another,  the  e plane  as  in  Figure  16,  the  points  correspond 
to  a mapping  in  which  points  of  a ray  from  the  origin  at  angle  9 from  the  real  axis  in 
the  e plane.  This  mapping  i,  not  conformal.  From  Equation  343  it  may  be  noted  that 
R ■ 0 at  the  origin  r * 0.  On  the  circumference  of  the  circles  R ■ 1 when  r * 1 . R 
is  complex  outside  the  circle  r > I,  Thus,  a Mach  circle  may  be  referred  to 
interchangeably  in  the  physical  or  t planes. 


i 

I 

I 


! 
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In  Equation  3*40  the  general  real-valued  solution  of  the  linearized  equation  was 
found  to  be  the  form  f(e)  + ft?).  Since  the  perturbation  velocity  components  u,  v,  and 
w are  solutions  of  this  type,  it  foltowi  that  they  are  the  real  parts  of  functions  of  the 
complex  variable  e and  complex  velocity  components  may  be  written 


1 u + iu* 

1 v + iv* 

1 w + iw* 


(344) 


where  u*.  v*,  and  w*  are  alio  real.  Since 


u*  - lu  - - \U(t) 
v*-  lv-  iK(«) 
w*  - iw  ■-  |W(«) 

The  components  u*,  v*,  and  w*  are  perturbation  velocity  components  of  some  other 
conical  flow. 

The  complex  velocity  components  are  related  by  the  compatibility  relationships 


If  V W' 

2Ttan*i“_(1  + e2)  " - i(i  - e2) 


(345) 


The  derivation  of  the  compatibility  relationships  is  included  in  Appendix  A. 

To  determine  the  flowfield  about  an  airfoil,  some  terms  and  nomenclature  must 
first  be  defined.  An  edge  may  be  leading,  trailing,  or  axial  which  Milne-Thomson 
denotes  as  I.,  T,  and  A.  If  a line  drawn  parallel  to  the  freestream  velocity  vector 
enters  the  airfoil  when  it  crosses  an  edge,  the  edge  is  culhd  a leading  edge.  If  the  line 
leaves  the  airfoil  as  it  crosses  the  edge,  the  edge  is  trailing.  If  the  edge  is  parallel  to 
the  freestream  velocity  vector,  it  is  axiul.  In  addition,  an  edge  may  be  termed 
subsonic,  sonic,  or  supersonic,  depending  upon  the  normal  component  of  the 
freestream  velocity,  A subsonic  edge  will  fall  within  the  Mach  cone  while  a supersonic 
edge  will  fall  outside.  A sonic  edge  corresponds  with  tl.e  Mach  cone.  Milne-Thomson 
denotes  this  by  subscripts  i,  s,  and  o for  inside,  sonic,  and  outside.  Thus,  for  example, 
an  intersection  may  be  denoted  L0T,,  indicating  a leading  edge  outside  the  Mach  cone 
and  a trailing  edge  inside.  For  an  axial  edge,  the  subscript  is  dropped  since  it  must  be 
subsonic. 

To  determine  boundary  conditions,  the  physical  situation  is  reduced  to  a Mach 
circle  in  the  physical  plane.  Consider  a LQL  airfoil  oriented  at  a small  angle  of 
attack,  a,  to  freestream  velocity,  as  shown  in  Figure  17a.  The  edges  8.  and  8,  are  at 
angles  5,  and  6,,  respectively,  to  the  axis  and  intersect  an  arbitrary  plane,  P, 
perpendicular  to  tne  axis  at  points  Lj  and  L2.  The  intersections  with  the  plane  P are 
shown  in  Figure  17b  with  corresponding  identification.  With  supersonic  edges,  two 


i 


(b) 


FIGURE  11.  (•)  PfcyriMl  iUyMwnUdun  of  L L Airfoil;  (b)  L L Airfoil 
oil  Pby#c*l  Hat  P. 


NWC  TP  S9S 1 


planes  may  be  drawn  from  the  edges  tangent  to  the  cone  at  points  A and  A'.  The 
points  at  which  the  Mach  cone  intersects  the  plane  of  the  airfoil  are  denoted  D(  and 
D,.  The  angles  between  the  tangency  points  and  the  plane  of  the  airfoil  are  denoted 
by  X. 


cos  xi " xc; " '<5T  at;  “ tan  cot  fil 


(3-46) 


Similarly,  cos  fi2  * tan  /u  cot  fi2. 

The  physical  significance  of  the  tangent  planes  may  be  realized  by  considering 
each  point  on  the  leading  edge  as  a generator  of  a Mach  cone  as  shown  in  Figure  18, 


MACH  CONU 


FIGURE  18.  Definition  of  Region  Where  TwoJJtmenienel  Bounded  Condition!  Apply. 


f 


l 


I 
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The  planes  formed  by  these  cones  are  the  planes  A|L,  and  A',L.,  Thus,  these  planes 
define  the  region  of  influence  of  the  airfoil  where  the  perturbation  velocities  are 
defined  by  two-dimensional  flow,  Therefore,  the  streamwise  components  of  the 
perturbation  velocities  in  Figure  17  are: 

u ■ u,  on  arc  D,A, 
u * - u(  on  arc  D(  A', 
u * u2  on  arc  D2A2 
u »-  u2  on  arc  D2A2 
u * 0 on  arc  AjA2 
u ■ 0 on  arc  A'j  A2 

where  u,  and  u2  are  appropriate  two-dimensional  perturbation  velocities 

The  case  of  an  L Tj  airfoil  is  shown  in  Figure  19a.  Since  the  leading  edge  does 
not  intersect  the  plane  of  the  Mach  circle,  the  ray  Is  extended  to  intersect  at  point  L'. 
Also,  since  the  trailing  edge  is  inside  the  Mach  cone,  the  tangent  planes  are  not 
present.  This  case  results  in  the  Mach  circle  representation  of  Figure  19b.  The 
perturbation  velocities  are: 

u ■ u0  on  arc  D2A 
u ■-  u0  on  arc  D2A'  ' 
u “ 0 on  arc  A'Dj  A 

These  cases  illustrate  the  essential  representation  for  other  cuses. 

The  basic  problem  in  conical  flow  is  to  find  a function,  fie),  analytic  within 
the  Mach  circle,  y,  whose  real  part  takes  on  specified  values  on  the  circumference. 
Since  a point  on  the  circumference  can  be  specified  by  the  angle  0 from  the  real  axis, 
the  real  purt  of  fie)  must  take  on  the  value  F(0)  on  the  circumference  where  F(0)  is 
given,  For  u point  a on  the  circumference 

a - a10  (3-47) 

or 

0 - (log  o)/i  (3-48) 


Thus, 


F(0)  - F((!og  o)/i|  ■ g(o)  (3-49) 

where  the  function  g(o)  is  likewise  known.  The  real  part  of  fie)  is  l/2fie)  + l/2f(e). 


I 


i 


i 
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On  the  circumference 


e = o and  a = e"ie  = 1 la 


therefore 


f(»)  + T(l/o)  “ 2g(a) 


(3-50) 


Multiplying  this  equation  by  do/2iri(o  - e)  and  integrating  around  the  Mach  circle,  y 


1 1 f ai/2242.  1 ftmi 

2m  Jy  a - e 2m  a - e m o - e 


(3-51) 


By  Cauchy’s  formula,  the  first  integral  is  f(e).  The  second  integral  may  be 
evaluated  by  the  residue  theorem.  Since  f(e)  has  been  assumed  analytic  within  the 
Mach  circle,  it  may  be  expanded  in  an  infinite  series. 


f(e)  ■ a0  + aj€  + a2e2  + • • • 


T(  1 /o) « ap  + a,  /a  + a2/o2  + 


(3-52) 


By  Cauchy’s  residue  theorem 


1 /*  5 Odo 


(3-53) 


-f 

Vy  on(o- 


(o-  e) 


1 0 for  n > 0 


(3-54) 


That  is,  the  residue  is  the  appropriate  term  of  the  Laurent's  series.  Thus,  from 
Equation  3-51 


+ 7(0)oit  f S&da 
' ' in  Jy  a - e 


(3-55) 


1 


ir 


i 


< 
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In  many  applications,  g (a)  is  constant  over  arcs  of  the  circumference.  Consider 
for  example  that  g (a)  = uQ  from  6 = - X to  + X and  is  zero  on  the  rest  of  the 
circumference.  Then, 


o - e 


(log  (0  - e)] 


0**! 


(3-56) 


for  e,  * and  * e“*\ 


Therefore, 


f(e)  - ■■  7(0)  + ^[log  (e,  - e)  - log  (e,  - e)l 


(3-57) 


J 


I 

I 

I 


i 

j 


7(0)  is  a constant  determined  by  other  information.  For  f(e)  ■ u0  when  e * 1 
or  0 ■ 0,  then 


u0  a - 7(0)  + t^-[log (e,  - 1)  - log (e,  - 1)] 


which  gives'" 


7(0)  - u0X/jr 


I 

The  same  principle  car  be  used  to  determine  a function  f(e)  analytic  within  a 
semicircle  of  unit  radius  where  ffe)  is  real-valued  on  the  diameter  and  the  real  part  of 
fi(e)  takes  given  values  on  the  semicircle.  To  solve  this  problem,  the  circle  is  complete 
by  reflection  as  shown  in  Figure  20.  The  function  f(e)  and  Its  reflection  are 
determined  to  be  complex  conjugates,  This  reduces  the  problem  to  the  one  Just  solved. 

! 


*!  * Thli  mult  li  baud  on  the  convention  that  the  branch  cut  !i  on  the  negative  real  axis  and  that  all 

| anjfes  thus  may  take  on  value  betwoon  ~v  and  <r. 

j 39 
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PIGUKE  20.  He  faction  of  Sm+CIrcfe. 


To  investigate  the  flowfleld  of  a two-dimensional  inlet,  the  model  of  a 
rectangular  airfoil  may  be  employed.  This  is  equivalent  to  a L A airfoil  in 
Milne-Thomson’s  notation.  The  physical  situation  is  shown  in  Figure  21a,  while  the 
projection  in  the  e plane  is  shown  in  Figure  21b.  The  boundary  conditions  are  also 
included  in  the  figure.  On  arc  DA,  the  perturbation  velocity,  u,  is  equal  to  the 
two-dimensional  perturbation  velocity  uQ  B oV  tan  p cosec  X.  For  5 * jt/2,  u0  is  equal 
to  aV  tan  p.  On  the  lower  arc  DA',  U is  -u0.  On  the  surface  of  the  airfoil  AD,  the 
imaginary  component  u*  of  the  complex  perturbation  velocity,  U , is  zero.  Consider  an 
airfoil  at  angle  of  attack,  a.  For  small  angles  where  sin  a * a and  coa  o ■ 1,  the 
freestream  velocity  may  be  considered  as  having  velocity  components  V,  0,  and  aV  in 
the  coordinates  of  Figure  14.  On  the  surface,  since  the  air  cannot  flow  through  the 
airfoil,  w + aV  ■ 0.  From  the  compatibility  relationship  of  Equation  3-45, 


If  W^_ 

3etanp  " _ 62^ 


II 

l' 


i 

l 

l 


I 


( 

I 

\ 

( 
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On  tho  surface  of  the  airfoil,  e is  teal  and  hence  (1-  e2)/2e  tan  /i  is  real.  Collecting 
the  real  terms  of  Equation  3-58, 

Since  w ■ -aV  * constant,  dw  ■ 0.  Therefore,  du*  * 0 or  u*  ■ constant. 
Without  loss  of  generality,  the  constant  may  be  considered  to  be  zero  so  that 


u*  ■ 0 


and 


u - U(e) 
on  the  airfoil  surface. 

To  solve  the  problem  of  the  L0A  airfoil,  define  m2  ■ e mapping  the  Mach 
circle,  7,  on  a semi-circle  in  the  v plane  which  is  completed  by  reflection  as  shown  in 
Figure  22.  If  a is  a point  on  the  circumference,  7' 

U(o)  + £/(  1 la) " 2u(o)  (3-59) 


noimi  n.  l0a  Mm  * nm. 
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As  shown  previously,  this  leads  to 


Since  u ■ u*  ■ 0 at  u ■ 0,  (/( 0)  ■ 0,  and  hence  V( 0)  ■ 0,  which  yields: 


iu0  (v,  + v)(t>.  + u) 


(3-61) 


Since  Milne-Thomson  is  concerned  with  the  pressure  forces  on  the  airfoil,  he 
does  not  determine  the  other  perturbation  velocity  components. 


EXTENSION  OF  THEORY 


The  compatibility  relationships  may  be  employed  to  determine  these  other 
components.  Equation  3-61  may  be  rewritten 


U(v)m-  -Ji  (log  (Dj  +V)  + log(P,  +!>)-  log  (V , - v)  - log  (P,  - M)]  (3-62) 


The  compatibility  relations  are: 


JML 

“twu  (1+e2)  i(l-e2) 


In  terms  of  the  variable  v where  v1  ■ e 


wf1-?1!-*/! 


NWC  Tl»  S95I 


Similarly 


V(v)'  • V(vY  gj 
W(eY  ' W(i>)'  jjj 

Thus  the  compatibility  relationships  become 

dg.....JfeX  to 

2e  tan /4  de  (j +e2^de  i(]-«2)de 
or  dividing  by  g- 

m u . m' 

tan?  (i+«2)  1(1  ~e2) 


In  terms  of  v 


Thus, 


and 


If 


W' 


2v2  tan m 0+»'4)  i(l  - »A) 


y . _ y 

2v2  tan  n 


JO  if 

2v i tan  n 


From  Equation  3-62 


iUfl  [ 1 « 1 t 

" ■ - T I v-7 + jpri  * + 'r' 


Using  the  compatibility  Equation  3-65, 

SjjHVU  I j — 4 — ! — + _) — 

» ‘T’-'’  •’i-"  T-' 

^ ■ irrars  (?♦  ‘,!)[p7T5 * ispn; + :7b + f~7 


(3-63) 


(3-64) 


(MS) 


(3-66) 


(3-67) 


(3-68) 
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Integrating  and  simplifying  yields: 

iu0  r 2 2 i /pi+p\  i /p . 

- ^log(”)~  j^log(  >o«(v,  +») 

+ P^  log  (P,  + p)  - v\  log  (v,  v)  - Pf  Ion  (P,  - v)j  + C (3-69) 


where  C is  the  constant  of  integration. 


To  evaluate  the  constant  of  integration,  the  boundary  condition  v » 0 at 
vm  1 +0i  ■ 1 is  employed.  Recognising  that  P,  ■ 1/V,  and  »»,  ■ 1/P(  and  evaluating  at 
v “ 1 


or 


Vm 


lu0 

3tr  tan  n 


-4P 


t 


+ P?  [log C^|  + 1)-  log(v,  - l)  + log(P,  +1)-  log (P , “ 1)] 

+ (log  (P,  + 1)-  log  (P,  •-  l)  + 'og(v,  + 1)  log (i/  - 1)]|+C 


V ■ 


IUq  ( _ , O',  + D(P,  + 1) 

5rTSHMp,'i + pi> + (v\ + p?> log v;-  n(v-  n 


+ C 


In  the  general  case, 

Vj  *0  + 61 


(3-70) 


and 


P,  -a~6i 

where  a ond  b are  real  constants. 

Substituting  the  general  values  for  v,  and  P,  and  recognizing  that  sinoe  v,  is 
on  the  circle  y\  ir  + b2  • 1,  Equation  3-70  becomes: 


iuo 

ymiwiart 


-8a  + (2flJ-262)logfrf 


+ C 


\ 
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The  phytical  perturbation  velocity  v is  the  real  part  of  the  complex  perturba- 
tion velocity  V, 

For  0<a<  1,  the  semicircle  In  the  v plane  which  repreients  the  circle  y in  the 
< plane  and  is  the  physical  part  of  the  problem,  the  term  1 + a/1  - a is  always  real 
and  positive.  The  singularity  at  a ■ 1 corresponds  to  a physical  situation  where  the 
angle  between  the  extension  of  the  leading  edge  of  the  airfoil  and  the  axial  edge  is 
equal  to  the  Mach  angle.  This  is  the  limit  where  the  model  of  an  L0A  airfoil  may  be 
applied.  Since  the  flwction  1 + a/1  - a is  always  real  and  positive  in  the  region  of 
interest,  the  angle  that  the  function  makes  with  the  real  axis,  which  it  the  principal 
value  of  the  imaginary  component  of  the  logarithm,  is  aero.  Thus,  the  only  real  part 
of  the  complex  perturbation  velocity  occurs  in  the  constant  of  integration.  Since  v ■ 
0,  at  v ■ 1,  the  real  part  of  the  constant  it  zero  and  the  imaginary  part  may  be 
ignored  since  it  will  not  be  required  in  determining  the  phytical  perturbation  velocities. 

Hence,  for  the  purpose  of  evaluating  the  physical  perturbation  velocity, 


V 


_ luo 
2?  tan  H 


> i 


i 


■i 


2a»,sr  - 2P,i'  + v]  log  O',  ■«■>») 


+ Pf  log  (P,  + V)  - >»!  log  (v,  - V)  - p]  log  (P|  - v) 


(3-71) 


From  the  compatibility  relationships,  as  shown  In  Equation  3-66: 


i 


b 


S 


W'  ■ - 
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Integrating  as  before 


jj2  ,0*  (V")'  ji  ‘°C  (V")  - £1 ,0*  (V“)-  *1  l0*  <V1  + •»> 

log  (P,  + v)  + log  (v,  - v)  + pj  Jog  (P,  - »oJ  + C (: 


(3-73) 


where  C it  the  constant  of  Integration. 

For  the  boundary  condition  at  v ■ l,wB0 

^""STarTillp!  ,0*^l  + + lo*(pi  + 1>“  ?|  " D 

“ vf  log  (P,  - 1)  - v\  log  (Mj  + 1 ) - P^  log  (P|  + 1)  + v\  log  (p,  - 1 ) 

+ P?log(P,-  1)1  +C 


Uq  ( , , (P(  " 1)(P.  + 1)1 

w " -2?Sui r i ' p? ) ,og ^ + w! -i)  i+c 


for  t>|  ■ a + hi  and  P,  ■ a - b i 


H,--2?T^!4fl6,ll0*(-1))i+c 


2nbUa 

H,-t^+C 

The  physical  perturbjtion  velocity  w is  the  real  part  of  the  complex 
perturbation  velocity  W 

'iabyXf. 

W ■ — " 1 ■ 4.  C* 

tan  p c 

where  C'  is  the  real  part  of  constant  C. 
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Since  w * 0 at  p * 1 


C' 


2a6u0 

taiuT 


Hence 


w"“3ra7i[-v“i^+2^v+2P^+^l08(^) 

- ^ log  (V")  ' pf  1°8  (’ V")"  A + *>  “ P1 

+ p?  log  (p,  - V)  + p?  log  (P,  - p)|-2  55~ 


+ =Ilog('iT") 


log(Pl  + p) 


(3-74) 


These  expressions  for  perturbation  velocities  have  been  evaluated  around  and 
within  the  circle  y'  for  an  arbitrary  value  of  p(  “0.8  + 0.61  which  corresponds  to 
5 ■ 116  degrees  ut  M ■ 2,  The  physical  perturbation  velocities  on  the  circle  y'  are 
presented  in  Figure  23.  The  magnitudes  and  signs  of  the  perturbation  velocities  are 
determined  by  the  equations  without  requiring  separate  determinations  of  the 
two-dimensional  perturbation  velocities  v0  and  wQ.  Consider  the  physical  situation  of 
the  thin  plate  at  positive  angle  of  attack.  The  two-dimensional  perturbation  in  the 
direction  of  the  flow  would  be  negative  on  the  lower  surface  as  the  velocity  is  reduced 
from  the  freestream . value.  On  the  upper  surface,  an  expansion  results  In  a higher 
velocity  or  positive  perturbation.  The  vortical  velocity  perturbation,  w,  will  have  the 
same  magnitude  and  direction  on  both  upper  and  lower  surfaces  and  will  be  negative 
or  downward  fur  positive  angle  of  attack.  The  sideward  velocity,  v,  may  be  expected 
to  change  direction  on  the  upper  and  lower  surfaces  and  to  be  positive  or  outward 
with  swoepback  toward  the  axial  edge.  For  the  case  of  no  sweepback,  6 * 90  degrees, 
the  two-dimensional  perturbation  v0  is  zero.  Thus,  the  perturbation  velocity  v is  zero 
everywhere  on  the  circle  y',  However,  due  to  the  imaginary  component,  v has  non-zero 
values  within  the  circle. 

The  perturhutlon  velocities  v and  w each  have  a singularity  at  the  point  P“0 
which  correaponds  to  the  axial  edge.  This  is  due  to  the  tenn  2/s» t p +•  2/Vj  m whose  real 
and  imaginary  components  approach  “ us  p approachea  zero,  There  appears  to  be  an 
additional  contribution  from  terms  of  the  form  log  [(Vj  * p As  p approaches  zero, 
the  real  part  of  the  logarithm,  which  is  equal  to  the  tog  of  the  amplitude,  approaches 
oo.  However,  examination  of  the  equations  for  perturbation  velocities  reveals  that  eaoh 
of  these  terms  may  be  combined  with  a similar  term  to  be  of  the  form  log 
((p,  + p)/(P|  - p)].  This  term  approaches  log  (1)  as  p approaches  zero,  The  term 
2/p,p  + 2/P(p,  however,  still  causes  a singularity  at  p « 0. 
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The  physical  perturbation  velocities  are  shown  in  Figure  24  as  functions  of 
location  in  the  v plane  for  the  arbitrary  value  of  « 0,8  + 0.6i.  The  singularities  in  v 
and  w are  apparent  for  small  values  of  radius, 

There  arc  also  singularities  when  v equals  v,  or  i», . At  this  point,  the  terms  in 
perturbation  velocities  v and  w of  the  forms  log  [(y,  - v)/v)  or  log(v,  - v)  become 
log(0)  which  cannot  be  evaluated.  In  addition,  the  perturbation  velocity  u has  a 
singularity  at  v,  and  P|  due  to  the  term  log  (s», +»»)(?,+»»)/(»»,  ••  v)(P,  - v)  which 
approaches  log  (<*), 

Finally,  singularities  exist  in  perturbation  velocity  u on  the  circle  7'.  The 
velocity  u is  proportional  to  the  imaginary  component  of  a logarithm.  The  principal 
value  is  the  angle  of  the  argument  of  the  logarithm.  The  real  part  of  the  argument  is 
negative  while  the  imaginary  part  it  zero,  Due  to  round-off  error,  the  imaginary  part 
has  a very  small  but  non-zero  value  which  may  be  either  positive  or  negative.  Thus, 
the  angle  is  evaluated  ns  plus  or  minus  tr  radians.  This  singularity  causes  no  practical 
problem  in  evaluation  of  velocities  since  this  will  occur  only  when  the  point  In 
question  is  exactly  on  the  circle  7'. 


Chapter  4.  STREAMTUBE  CALCULATIONS 

To  evaluate  the  amount  of  sidespill  and  the  additive  drag  of  a two-dimensional 
inlet,  a trace  of  the  captured  atieamtubc  was  employed.  The  velocities  within  the 
conical  flow  region  cun  be  calculated  by  Equations  3-62,  3-71,  and  3-74,  Within  the 
two-dimensional  flow  region,  the  perturbation  velocities  are  also  known,  To  determine 
tne  captured  stream  rube,  a numerical  calculation  of  streamlines  from  the  inlet  cowl 
was  performed.  A computer  program  wbr  written  in  FORTRAN  V for  the  Univac 
I ilO  computer.  This  program  is  listed  In  Appendix  C and  a flow  chart  is  presented  in 
Figure  25.  The  numerical  computer  soiutton  was  selected  for  several  reasons.  The 
repetitive  nature  of  the  calculations  In  tracing  the  streamlines  from  the  cowl  to 
freestream  permits  the  application  of  a computer  to  this  problem.  This  permits  the  use 
of  small  incremental  distunces  to  accurately  portray  the  streamtube.  Determination  of 
additive  drag  is  also  performed  on  the  computer  by  numerically  integrating  the 
pressures  over  the  surface  of  the  streamtube,  in  addition,  the  computer  simplifies  the 
solution  of  the  complicated  and  complex  equations  of  the  conical  perturbation 
velocities.  Once  a variable  has  been  identified  as  a complex  variable,  the  equations  In 
FORTRAN  are  written  as  with  real  variables,  Arithmetic  equations  are  performed  using 
the  complex  variables,  Functions  such  as  square  roots  and  logarithms  must  be 
identified  as  complex  functions.  Hand  calculation  of  perturbation  velocities  are  tedious 
und  prone  to  error. 
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FIGURE  25.  (Contd,). 


DESCRIPTION  OF  COMPUTER  PROGRAM 


To  determine  the  perturbation  velocities  at  any  point,  the  inlet  ramp  is 
represented  as  a thin  airfoil  in  the  xy  plane  at  u negative  angle  of  attack  a.  This  is 
accomplished  by  a rotation  transformation  of  the  ramp  and  cowl  coordinates.  The 
rotation  changes  the  freestream  velocity  components  to  V,  0,  -oV.  I he 
two-dimensional  perturbation  components  are: 

u0  =*  aV  tan  m cosec  X (4-1) 
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where  cos  X * tan  n cot  5 

vQ  = oV  tan  n cot  5 cosec  X 
w0  » aV 

Derivation  of  these  velocity  components  is  included  in  Appendix  D. 

Both  the  conical  and  two-dimensional  shocks  may  be  represented  as  Mach 
waves  in  the  rotated  coordinate  system.  For  small  angles  and  Mach  numbers,  the 
two-dimensional  shock  angle  is  approximately  equal  to  the  Mach  angle  plus  the 
deflection  angle  as  illustrated  in  Figure  26,  The  conical  shock  wave  angle  is 
approximately  equal  to  the  Mach  wave  angle  measured  from  the  axis,  over  the  same 
range  of  Mach  number  and  conical  deflections  as  shown  in  Figure  27. 

To  trace  the  streamtube,  initial  points  are  selected  uround  the  cowl  lip  and 
sideplate  leading  edge.  At  the  side  edge  of  the  ramp,  there  are  singularities  in  both  v 
and  w.  To  avoid  this,  two  points  were  chosen  in  the  neighborhood  of  the  singularity, 
one  on  the  vertical  sideplate  and  one  on  the  ramp  surface  slightly  within  the  inlet. 
From  the  latter  point  the  vertical  velocity  is  zero  and  the  streamline  is  constrained  to 
the  ramp  surface.  Thus,  this  defines  the  lower  edge  of  the  streamtube.  The  presence  of 
the  singularity  also  makes  it  necessary  to  employ  small  differences  in  axial  location  for 
the  streamline  trace  in  this  region.  A large  difference  results  in  a vertical  change  such 
that  the  streamline  is  projected  through  the  ramp. 

In  o^der  to  calculate  perturbation  velocities,  a determination  must  first  be  made 
of  the  type  of  flow  existing  at  any  point.  Three  types  of  flow  were  considered: 
conical  flow,  two-dimensional  flow,  and  undisturbed  flow.  I'he  square  of  the  physical 
radius  y2  + z2  is  compared  with  the  square  of  the  rudius  of  Mach  cone,  x2  tan2  /u,  to 
determine  if  a given  point  is  within  the  cone  and  hence  in  the  region  of  conical  flow. 


If  the  point  is  in  this  region,  the  position  is  recalculated  in  the  form  of  Buscmunn's 
complex  variable,  e,  and  then  by  the  complex  square  root,  as  »>.  The  conical 
perturbation  components  are  then  evaluated. 

If  the  point  in  question  is  outside  the  cone,  the  angle  tan" 1 z/x  is  compared 
with  the  Mach  angle,  n,  If  the  angle  is  less  than  n,  the  point  is  within  the  region  of 
two-dimensional  flow,  and  the  perturbation  velocities  are  two-dimensional  components. 
If  the  angle  is  greater  than  n,  the  perturbation  components  are  zero. 

Once  the  perturbation  velocity  components  have  been  determined,  they  are 
added  to  the  freestream  velocity  yielding  the  local  velocity  at  each  point.  The  velocity 
vector  is  projected  an  incremental  distance  upstream  to  a new  set  of  spatial 
coordinates.  A rotational  transformation  converts  the  streamline  coordinates  hack  to 
the  physical  coordinate  system  for  printing.  This  procedure  is  repeated  for  each 
streamline  and  then  for  a new  series  of  incremental  distances  until  the  streamtube  is 
determined  back  to  the  inlet  ramp  leading  edge.  The  order  of  calculation  could  have 
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FIGURE  27,  Approximation  for  Conical  Shock  Wavei. 
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been  different  such  that  each  streamline  wus  traced  to  the  leading  edge  before  a new 
streamline  was  examined.  The  procedure  chosen  yields  data  groupings  of  cross  sections 
at  approximately  constant  axial  stations, 

The  static  pressures  ore  evaluated  at  each  point.  From  the  perturbation  velocity 
components,  the  local  velocity  components  are  determined.  For  simplicity,  all  velocities 
have  been  ratioed  to  the  freestream  velocity,  V.  Initially,  the  ratio  V/a*,  where  a*  is 
thw  speed  of  sound  at  sonic  conditions,  has  been  determined  from  the  freestream  Mach 
number. 


V/a*„|i^iM2(i  M2r‘ 1/2 


(4-2) 


Since  the  compression  process  is  adiabatic,  a*  remains  constant.  Thus,  when  the 
magnitude  of  the  local  velocity,  Vg,  is  determined,  the  ratio  Vg/a*  may  also  be 
determined. 


Vg/fl*-(Vg/V)X  (V/a*) 


(4-3) 


From  the  value  of  Vg/a*,  the  ratio  of  local  static  to  total  pressure  may  be 
determined* 


(p/pt)B  “ [l  ' ] (4-4) 

For  the  calculations  herein,  it  was  assumed  that  ptg  ■ pt-1  which  is  reasonable 
for  the  small  ramp  angles  and  low  Mach  numbers  considered.  It  would  be  simple  to 
include  the  ratio  of  total  pressures  across  the  shock  wave. 

Once  the  streamline  has  been  traced  upstream  to  the  ramp  leading  edge,  the 
final  cross-sectional  area  is  evaluated  by  numerical  integration  and  divided  by  the 
capture  area  to  yield  the  mass  flow  area,  A0/Ao,  for  both  spillage  and  below  design 
Mach  number  operation. 

Following  this,  a numerical  integration  of  the  static  pressures  over  the 
streamtube  is  performed  to  obtain  the  additive  drag.  The  static  pre  sure  has  been 
calculated  at  each  printed  point  of  the  streamtube.  A numerically-averaged  pressure 
over  the  projected  area  of  each  element  is  summed  to  yield  the  drag  force  which 
determines  the  additive  drag  coefficient. 
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* Rtfti  to  Equation  S3  In  Equation!,  Tahiti,  and  Cham  for  Compnathk  Flow,  by  Am#i  Rimrch  Staff. 
National  Advtiory  Commit*#*  for  Aeronautici,  Moffatt  Field,  Calif.,  1953.  (NACA  R#port  1135.) 
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„ (P-P0WA  £C(P/Pt)«  " P0/PtJ  AA 

Cd.‘  9Aref  - (4*5) 

Thus,  this  computer  program  hat  yielded  both  the  mau  ipltlage  and  the 
additive  drag  coefficient  from  the  email  perturbation  velocity  components. 

For  presentation  purposes,  the  computer  program  was  coupled  with  a plotting 
routine,  DISSPLA,  developed  by  Integrated  Software  Systems  Corporation.  Since 
plotting  routines  are  generally  unique  to  each  computer  system,  further  description  and 
listing  of  plotting  call  statements  are  omitted. 

Typical  results  of  the  program  am  included  in  Figures  28  through  33.  The  inlet 
considered  has  a 10-degree  ramp.  The  cowl  is  located  10  units  of  length  from  the 
leading  edge  with  a height  of  5 units,  inlet  width  is  10  units.  Freestream  Mach 
number  is  2.2.  The  first  five  plots,  Figures  28  through  32,  show  the  shape  of  the  cross 
section  of  the  streamtube,  from  the  mldpiane  to  the  axial  edge  at  various  longitudinal 
locations  from  the  cowl  lip  to  the  leading  edge.  Superimposed  on  these  plots  are  the 
boundaries  of  the  two-dimensional  and  conical  flow  regions  to  illustrate  the  analytical 
model  for  each  portion  of  the  streamtube.  At  the  leading  edge,  Figure  32,  the  flow  is 
undisturbed  and  there  are  no  regions  of  perturbation.  Integration  of  this  cross  section 
yields  the  mass  capture  area  ratio  of  the  inlet.  Figure  33  presents  a perspective  view  of 
the  inlet  and  streamtube. 


PAWLIKOWSKI’S  ANALYSIS 

At  this  point,  it  Is  worthwhile  to  examine  the  investigation  of  Pawlikowski 
(Ref.  10).  Pawlikowski  has  performed  an  analysis  of  the  perturbation  velocities  similar 
to  that  already  discussed.  To  determine  the  velocities,  he  employed  a technique 
developed  by  Lagerstrom  (Ref.  13).  This  method  will  be  described  here  in  the 
coordinate  system  of  Milne-Thomson  shown  in  Figure  14  rather  than  In  Lagerstrom’s 
notation.  Starting  from  the  small  perturbation  equation, 


3x2  dy2  9a2 


(4-6) 


a Mach  number  of  \/T  ia  assumed.  This  reduces  the  term  M2  - 1 to  unity.  Other  Mach 
numbers  are  treated  by  distorting  the  longitudinal  dimension  by  a Prandtl-Glauert 
transformation. 


The  imall  perturbation  equation  is  further  reduced  by  expressing  it  in 
cylindrical  polar  coordinates.  Conical  flow  is  evoked  by  observing  that  the  flowfleld  ia 
completely  apecifled  if  it  ia  known  within  the  Mach  circle  in  the  plane  x * 1 . In  this 
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LL 


plane,  with  a freestream  Mach  number  of  ^T,  the  Mach  circle  has  a radius  of  unity 
since  tan  ju  " I.  In  the  plane  x ■ 1,  the  small  perturbation  equation  may  be  written:* 


. i jl  as.  as  . 1 £s „ 


(4-7) 


This  may  be  reduced  to  Laplace’s  equation 

a2s  . i as  A i a^s 


3R2  + R + r2  “ 


by  a Tschapiygin  transformation 


1 -VI 


(4-8) 


which  is  identical  to  Equation  3-43.  As  in  Milne-Thomson’s  analysis,  Busemann'a 
complex  variable  e,  and  the  compatibility  relationships  for  M«v/T  or  tan/i*l  are 
introduced.  At  this  point,  Ugerstrom  determines  the  perturbation  velocities.  The 
boundary  conditions  for  a symmetrical  thin  body  in  the  xy  plane  are: 


w » w0  on  the  upper  surface 
w ■ -w0  on  the  lower  surface 
w * 0 In  the  xy  plane  off  the  body 


Consider  a general  edge  within  the  Mach  cone  intersecting  the  plane  x ■ I at 
e0  (similar  to  point  T in  Figure  19).  The  boundary  conditions  may  be  satisfied  by  the 
imaginary  part  of  n logarithmic  singularity  nt  e0.  Thus,  he  writes: 


iw„ 


W *-  — log(e-  e0)  + f(e) 


(4-9) 


where  Ke)  is  analytic  at  <0  and  purely  imaginary  on  the  real  axis.  From  the 
compatibility  relationships,  he  obtains 


t \u  w0/  2 W t \ . 2«  . 

W t (eJ  .if'o)  l(e2-  l/(e> 


(4-10) 


For  the  special  case  where  the  edge  is  parallel  to  the  flow  direction,  «0  * 0 and 
W has  a logarithmic  singularity  at  e0  “ 0 which  is  represented  by  -i(w0/n)  log  e, 


* Equation  1-38  In  Hafmnc*  13, 
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The  boundary  conditions  are  written  for  the  planar  body  with  axial  edge 
equivalent  to  the  L0A  airfoil.  These  boundary  conditions  are  shown  in  Figure  34.  The 
singularities  at  e,  and  e2  are  also  expressed  as  logarithmic  singularities.  Combining 
these  with  the  singularity  at  the  axial  edge  yields 


W*-—  [Jog(e-  ej)(<?-  e2)“  lo8«l 


(4-11) 


which  satisfies  the  boundary  conditions  of  Figure  34. 
Employing  the  compatibility  relationship 


(4-12) 


Placing  the  terms  over  a common  demoninator  and  simplifying  yields: 


\ J "‘"o 


(4-13) 


FIGURE  34.  Boundary  CanUHom  in  < Wane. 
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Thus, 


Since  e2  is  the  complex  conjugate  of  g)  on  a circle  of  radius  1,  g|g2  = 1. 


1 > 
£ 
o 

i 

n 

(e-  CjKe  - e2) 

J 

This  may  be  expanded  to  yield: 


or 


d V 2w0 


d U m 2wo 
HFa"T' 


(e- e2)-  (e- g]) 

(e  - e [ )(g  - e2!)(ei  - «2). 


(e-  - e2.)  + fiT-  e2)<!e2  - ef) 


(4-14) 


Since 


and 


e [ * e1?k  * cos  X + i sin  X 


*=  e-iX  « cos  X - i sin  X, 

il 

g(  - c2  “ 2i  sin  X 

Equation  4-14  may  be  rewritten  as 


dU  lwo 

de  trsln  X 

Integration  yields 

iwn 


1 


1 


<“<l  e-e2 


e - e-> 


^"nrirxlo8^7-  + c 


(4-15) 


(4-16) 


(4-17) 


The  constant  of  integration,  C,  is  evaluated  by  the  boundary  condition  at 


Wn 


u-R«(^)"uo*-3JiT 

By  geometry,  as  shown  in  Figure  35, 


Re 


"l“el 

■ll0?-nrr 


since  the  imaginary  part  of  the  logarithm  is  the  argument. 
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SUMMATION  op  ANGUS  Of  TWANOU  - ff 

lift  - )i*y 
9 ■ *X 

PIOURB  31.  G*om*trfc«i  SoIuMm  of  Constant  of 

Intapiboa. 


w0 

c“-ffjurxtff+xl 

w0  e - « 2 

rjm-*Wl  Hlogyr^  + ir  + M 


Differentiating  this  yields: 


dU 

iw0 

1 

nr"1 

ir  sin  X 

e-  «2 

This  is  combined  with  the  other  compatibility  relationship, 
dV  1 (r  + l\  df/ 

-avr  “ - *2  V“"T"j 


to  obtain 


dV 


*wo  e2  + 1 

-1 

2tr  sin  X e 

e - «i 

(4*18) 

(4-19) 

(4*20) 

(4-2') 
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liquation  4-21  is  rewritten  as 


dV  ‘w0 


-j — * * -*• 

de  if 


cos  X 1_ + cor.  X 1 i_ 

slnX  c - fi  sLi  X e - ej  e 


(4-22) 


This  equivalence  is  more  easily  verified  by  working  backwards. 
Recognizing  that,  for  the  complex  conjugates  ej  and  e2, 


+ e2  “ 2 cos  X 


and 


«fj  — e2  * 2i  sin  X 
Equation  4-22  may  be  written 


|^o  cosXf  _!_  i _ isinX 
de  **  " tr  sin  X ~ e - e ^ + e - e2  " e cos  \ 


dV  1*0 

(ei  +e2> 


1 • + 1 


c,  - 


e-  e]'  + e-  e2  “ etc , + e2) 


(4-23) 


Placing  this  expression  over  a common  denominator  and  simplifying  by 
subtracting  and  dividing  by  identical  terms  yields: 


dK 


iwr 


de  * “ 2i'sin  \ 


+<1<2  ~elg2  " *1*2 

e(e - e,)(e  - «,) 


dV 


lWn 


de"*  ~ 2tr  slnX 


r,2  i.e|e2)(e;-  et) 

1 e(e  - «jj)(e  - e2)  ^ 


(4-24) 


Since  e2  is  the  complex  conjugation  of  on  a circle  of  unity  radius, 


e,e2  - 1 


Hence 


dV 


iw0  e2  + 1 


37  * * 2v  sin  T e He  - e ; )(c  - «?  ) 


e2  ~ el 


lwo  " e2> 

2iir  slnX  e [ « - e,  )<e ■r*^f 
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iw 


dV  

37"-5FiST5:“7 


o e2  + l 


1 +•  1 


e - ej  e - c2 


which  is  identical  to  Equation  4-21. 

Since,  from  two-dimensional  linearized  theory  shown  in  Appendix  D, 
vQ  * -u*,  cot  5 cosec  X ■ Wq  tan  y cot  6 cosec  X 


and 


cos  X * tan  ju  cot  6 

. . cos  X 

v0  - w0  cos  X cosec  X ■ w0 


(4-25) 


From  this,  Equation  4-22  becomes 


dV  ivo 

le  — 7 

Hence,  by  integration 
iv, 


'e  - «,  + e - «2 


tre 


o • “ c2  wo 

Flo*rrj7-Tl0*etC 


(4-26) 


(4-27) 


The  constant  is  evaluated  as  before  with  the  boundary  condition  that  at  e “ - 1 , 
v ■ Re(  V)  - v0 


since 


and 


Re 


- * - e, 


Re  ( log  - 1 1 "log  1-11*0 


Therefore 


0-“(»r  + X] 


v0  e-c,  w0 

V m rf  [-  i log  + * + XJ  - log  e 


(4-28) 


(4-29) 


Thus,  the  complex  perturbation  velocities,  as  far  as  their  contributions  to  the 
physical  or  real  components  are  concerned,  have  been  determined  by  Lagerstrom  for  a 
freestream  Mach  number  of  \J~T.  These  components  are  further  reduced  to  real  velocity 
components  by  a geometrical  argument, 
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Consider  a general  point  e ■ Reifl  within  the  circle  of  unity  radius  as  shown  In 
Figure  36. 

Re(-i  log (« j - «))  "argCe,  - e)  — 0, 

Ret-i  log  (e2  - «)]  ■ arg  (e2  - «)  ■ $2  (4*30) 

Therefore  from  Equations  4-19,  4-29,  and  4-11 

u0 

u ■ Re(J/) m + ff  + X] 

V w 

v-  Re(F)  - ■"  [$2  - + log  R 

w«Re<KO-~l0-  (*,  + *,)]  (4-31) 

Pawlikowski  employs  these  perturbation  components  to  determine  the  amount 
of  sidespiil  and  the  additive  drag  coefficient.  He  considers  only  the  crosswise  velocity 
component  v, 

V w 

V--I02-0J +ff  + XJ-"“logR  (4-32) 


nOUM  34.  OMMMtrical  lUspaHOtM  of 

nrannH  twp**/ 
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Pawlikowski  first  considers  the  simple  cue  of  an  inlet  with  no  leading  edge 
sweepback.  For  this  caw,  the  two-dimensional  spanwise  velocity  component  v0  is  zero. 
Hence,  Equation  4-32  reduces  to: 


w„ 


v ■ - — log  R 

In  physical  coordinates  this  equation  may  be  written  as: 


(4-33) 


w„ 


V " - 


i-v/T7 


*~  log  - — *7 — (4-34) 

The  average  velocity  v between  any  two  radii  r2  and  r2  may  be  written  as 

vdr 


t 


ri 


From  this  equation,  Pawlikowskl  obtains: 


(4-35) 


w 

• Kr  '-T)  r log p==r+  sin*  ’ 

2 1 I 1 -\A-  r 


(4-36) 


For  the  case  where,  the  Mach  cone  intersects  the  cowl  of  the  inlet,  integration 
is  performed  from  r2  «0  to  r2  * 1.  For  this  case, 


w„ 


T)m 

4* 

To  determine  the  man  flow  rate  of  spillage,  Pawlikownki  assumes  that  the 
density  is  constant  and  may  be  evaluated  from  oblique  shock  relationships.  This 
assumption  is  questioned  in  view  of  the  wide  variation  in  crosswise  perturbation 
velocity  as  demonstrated,  for  example,  in  Figure  24.  Although  this-  figure  is  for  a ramp 
with  sweepback,  there  is  a singularity  at  r«0  for  any  value  cf  v, , including  Mj  ■ 
1 /&i  deg  which  represents  the  ramp  without  sweepback.  Thus,  a large  variation  in 
local  density  is  expected. 

With  the  auumption  of  constant  density,  the  maw  (low,  rit,  may  be  written  as: 


m-PvAlide 


(4-37) 


where  Alide  is  the  area  bounded  by  the  Mach  cone  and  the  cowl,  u defined  in  Figure 
37a. 
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In  keeping  with  standard  inlet  nomenclature,  the  spillage  mass  <low  may  be 
expressed  as  an  area  of  streamtube  at  freestream  condition 

m " PmVmA0SS 

where  A0gS  the  streamtube  area. 

Thus 


^oss.jLJ-  (4-38) 

P~^ 

From  gas  dynamic  relations  across  an  oblique  planar  shock, 


» ft  ♦ I) llni  8,  ,4.3,) 

P-  (7  - 1)  M2  sin2  fl,  + 2 

In  the  more  general  case,  where  the  Mach  cone  does  not  Intersect  the  cowl,  ai 
shown  in  Figure  37b,  the  area  of  spillage  Is  separated  Into  two  regions.  The  first 
region,  A,,  is  bounded  by  the  Mach  cone  and  is  equivalent  to  the  case  already 
described.  The  second  region,  A„,  is  bounded  by  a streamtube  parallel  to  the  ramp 
surface.  The  average  velocity  between  each  ray  may  be  determined  by  Equation  4-36 
and  the  average  over  the  area  A„  is  calculated  by  numerical  integration.  Thus,  the 
mass  spillage  for  this  cue  may  be  written  as: 


A0ss  p ^l^l 


*11*11 
^ re] 


(440) 


Altnough  Pawlikowski  uses  the  spillage  area.  A|idf,  as  a reference,  any  reference 
area  could  be  employed. 

This  portion  of  Pawlikowski’s  analysis  is  considered  to  be  in  error.  He  has 
made  an  a priori  assumption  concerning  the  shape  of  the  captured  streamtube.  He  has 
assumed  that  the  spillage  effect  is  a small  perturbation  to  the  two-dimensional 
analytical  model.  Thus,  the  density  is  assumed  constant  and  the  area  of  sidespill  is  the 
same  as  the  side  of  the  two-dimensional  streamtube,  The  analysis  already  presented 
herein  has  demonstrated  that  the  streamtube  resulting  from  spillage  has  a complicated, 
three-dimensional  shape.  It  is  difficult  to  identity  the  region  at  the  side  of  the  inlet 
through  which  spillage  occurs.  Thus,  the  assumption  of  the  two-dimensional  streamtube 
Is  considered  to  be  incorrect. 


Pawlikowski  then  evaluates  the 
Dt,  may  be  written  approximately  as 

Dt  ■ (Pi  - Pm)  Aoss 


additive  drag.  He  states  that  the  additive  drag, 

(441) 
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This  result  assumes  constant  static  pressure.  The  urea  AqSS  would  be  the 
cross-sections!  area  at.freestream  conditions  of  the  spilled  strcamtube  and  under  these 


assumptions, 


D,  -/( p - p*)  dA  • (ps  - p„)/dA  - (p,  - p„)A0 


(442) 


These  assumptions  are  questionable  for  the  same  reasons  previously  discussed 
when  a constant  density  was  assumed.  There  is  a large  variation  in  local  velocity  which 
implies  variation  in  local  physical  properties. 


Pawllkowski  expresses  the  additive  drag  in  coefficient  form: 
CD."qT* 


(443) 


where  A,  is  the  inlet  area. 


From  Equation  442 


Pj  _ Pm  A(j 


(444) 


The  first  term  p,  - p„/q  is  the  pressure  coefficient  behind  the  planar  shock 
wave.  Thus,  Pawllkowski  writes: 


Aoss  Aild* 


Cd«  c^..d.^r 


(445) 


Pawllkowski  continues  to  consider  more  general  inlets,  specifically  those  with 
leading  edge  sweepback  and  yawed  side  edges.  Tills  involves  more  difficult 
determinations  of  the  average  velocity.  However,  ail  of  his  analysis  assumes  constant 
properties  and  knowledge  of  the  stream  tube  shape  and  thus  is  considered  too 
simplified  for  correct  application. 


Chapter  5.  APPLICATIONS  OF  COMPUTER  PROGRAM 

The  computer  program  for  tracing  the  streumtubr  and  determining  the  additive 
drag  has  been  applied  to  several  cases  to  determine  the  effect  of  various  inlet 
configurations  and  conditions. 
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COMPARISON  WITH  EXPERIMENTAL  DATA 

This  computer  program  has  been  applied  to  the  analysis  of  an  inlet 
experimentally  investigated  by  Kemper  (Ref.  14).  This  inlet  is  illustrated  In  Figure  38. 
The  inlet  is  a two-ramp  design  with  an  initial  ramp  angle  of  9 degrees  and  a second 
ramp  angle  of  19.75  degrees.  A straight  sidewall  extends  from  the  cowl  lip  to  a point 
just  forward  of  the  second  ramp.  The  inlet  mass  flow  was  determined  for  this  case  and 
compared  with  measured  data,  To  model  the  inlet  analytically,  only  the  first  ramp  was 
considered.  Above  the  design  Mach  number  of  2,3,  the  planar  shock  from  the  second 
ramp  is  contained  wholly  within  the  sidewalls.  At  Mach  numbers  below  2.3,  the  region 
of  the  second  shock  forward  of  the  sidewall  is  small  relative  to  that  of  the  tint  shoot. 
Thus,  the  contribution  of  the  second  shock  is  assumed  to  be  small  and  has  been 
neglected.  The  uctuul  inlet  employs  both  sidewall  and  throat  bleed.  To  account  for  ti.e 
bleed  mass  flow,  an  estimate  of  the  bleed  was  determined  by  means  of  the  automated 
procedure  of  Ref.  8.  This  value  was  subtracted  from  the  analytical  determination  of 
streamtube  cross-sectional  area  to  yield  a net  mass  capture  urea  ratio  for  comparison  to 
measurements  with  both  spillage  and  bleed.  The  eompurlson  of  capture  area  ratio  is 
presented  In  Figure  39,  The  upper  curve  is  that  which  would  be  predicted  by 
two-dimensional  analysis  without  spillage  or  bleed.  Above  the  design  Much  number,  tie 
two-dimensional  capture  area  ratio  is  unity  or  fill!  capture,  Bleed  ar.d  sidespil!  tvdi:e 
this  value.  The  capture  urea  ratio  predicted  by  the  present  method,  modified  for  bletd, 
is  compared  with  measured  flow  rates  for  two  inlets.  The  two  inlets  hove  the  sane 
aerodynamic  shape  but  are  manufactured  by  different  processes.  Inlet  4 is  machined 
front  an  investment  casting  while  inlet  5 Is  formed  from  sheet  metal.  The  agreement  of 
the  present  method  with  the  data  is  reasonable.  The  capture  area  ratio  predicted  >y 
Powllkowski’s  method  is  ulno  presented  in  Figure  39.  The  value  of  spillage  in  addition 
to  the  bleed  flow  was  determined  by  the  automated  procedure  of  Ref.  8.  The  logic  of 
this  automated  procedure  docs  not  permit  calculation  of  spillugc  ut  or  above  the  design 
Mach  number  of  the  inlet.  Hence,  the  comparison  is  limited  und  no  conclusion  muy  be 
drawn.  No  determination  of  additive  drug  was  made  in  the  test,  amt,  hence,  10 
comparison  is  possible  for  the  predicted  value  of  udditive  drug. 


EFFECT  OF  INLET  DESIGN  VARIABLES 


The  effects  of  various  design  variables  were  examined  theoretically  by  he 
method  described  herein.  For  these  studies,  u simple  single  rump  two-dimensional  ir  let 
was  considered.  The  baseline  inlet  is  shown  in  Figure  40.  The  inlet  consists  o;  a 
10-degree  rump  which  is  10  units  in  length,  Since  in  both  conical  Dow  md 
two-dimensional  flow,  tho  How  depends  only  on  direction,  the  units  of  length  are  not 
essentia!,  only  the  ratios  of  lengths.  Thu  cowl  is  5 units  high  and  the  inlet  is  10  units 
in  width.  The  baseline  inlet  has  no  leading  od  ;e  sweepback  or  sidewalls, 
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Effect  ot  Sidcspillage 

The  results  of  the  three-dimensional  analysis  described  herein  were  compared  to 
those  of  a two-dimensional  analysis  based  on  linearized  theory.  The  effect  of 
sidespillaee  on  the  captured  streamtube  is  illustrated  in  Figure  41.  The  case  of  a 
10-degree  ramp  with  the  freestream  Mach  number  of  2.5  was  considered  for  this 
example.  The  streomtube  determined  from  the  three-dimensional  analysis  is  the 
irregular  area  ABCD,  The  upper  portion  of  the  streamtube,  line  AB,  reaults  from 
tracing  the  streamlines  upstream  from  the  top  surface  of  the  cowl.  Freestream  flow 


FIGURE  41.  Completion  at  Captured  Msaaintuhei,  M » 2.5,  • ■ 10  Degree*. 
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which  is  above  this  line  is  spilled  over  the  top  of  the  cowl  and  is  not  captured  by  the 
inlet.  The  side  portion,  line  BO,  is  determined  from  the  side  of  the  inlet.  Flow  to  the 
right  of  this  line  is  spilled  to  the  side  of  the  inlet.  Comer  B originates  at  the  cowl  of 
the  inlet  and  thus  divides  the  flow  over  the  comer  from  the  sidespillage.  The  rectangle 
AEFD  is  the  captured  streamtubc  for  the  twcKlimensional  model.  Note  that  since  the 
flow  is  two-dimensionrl  at  the  niidplane  of  the  inlet  for  this  case,  the  heights  of  the 
streamtubes  at  this  point  are  identical.  The  excess  spillage  over  the  two-dimensional 
case  is  the  difference  between  the  spillage  of  area  GEFC  and  the  gain  of  area  ABG.  It 
is  obvious  in  this  case  that  there  is  a loss  in  capture  area  due  to  sidespillage.  The 
difference  in  capture  area  ratio  is  presented  for  the  various  Mach  numbers  in  Figure  42 
for  a ramp  angle  of  10  degrees.  In  this  case,  the  net  sidespillage  is  an  almost  constant 
10%  over  the  Mach  number  range  considered.  The  net  sidespillage  will  vary  depending 
upon  the  Mach  number  and  inlet  geometry. 

The  effect  of  sidespillage  on  the  additive  drag  coefficient  for  the  same  case  is 
shown  in  Figure  4i.  The  effect  of  sidespillage  as  revealed  by  the  three-dimensional 
analysis  is  a significant  decrease  in  inlet  additive  drag  as  compared  with  that  predicted 
by  two-dimensional  theory.  This  results  from  a combination  of  three  effects. 
Examination  of  the  pressure  forces  has  revealed  that  the  upper  surface  of  the 
streamtubc  is  the  major  contributor  to  the  total  drag  force.  The  presence  of 
sidespillage  reduces  the  surface  area  on  the  upper  surface  since  some  of  the  flow  is 
spilled.  In  addition,  the  deflections  on  the  upper  surface  of  the  streamtube  are  less  for 
conical  flow  as  evidenced  by  the  upper  region  ABG  in  Figure  41.  This  results  in 
reduced  local  pressures  and  hence,  reduced  drag.  Finally,  the  presence  of  high  velocities 
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near  the  side  edge  of  the  ranip  can  result  in  negative  yiessure  coefficients.  In  many 
cases,  the  side  of  the  streamtube  1g  found  to  reduce  the  total  drag.  The  net  effect  of 
these  three  factors  is  a decrease  in  additive  drag  resulting  from  sidespiilage  as  compared 
with  the  two-dimensional  model.  This  result  causes  doubts  about  the  method  of 
Pawlikowski  in  which  it  is  assumed  a priori  that  the  spillage  drag  is  an  additional  term 
to  the  two-dimensional  additive  drag. 


Effect  of  Ramp  Angle 


The  first  design  variable  examined  was  the  angle  of  the  inlet  ramp.  Deflections 
of  5,  10  and  15  degrees  were  considered  for  Mach  numbers  of  2.0,  2.5,  and  3.0.  The 
effect  of  ramp  angle  on  the  inlet  capture  area  ratio  is  presented  in  Figure  44.  As 
expected,  the  capture  area  ratio  decreases  with  increasing  ramp  angle  for  all  Mach 
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numbers,  The  freestream  shapes  of  the  captured  streamtubes  are  shown  in  Figure  45 
for  a freestream  Mach  number  of  2,5.  Similar  trends  are  observed  at  the  other  Mach 
numbers,  As  may  be  noted  from  the  figure,  both  the  spillage  over  the  cowl  and  the 
sidespillage  are  increased  with  ramp  angle,  decreasing  the  capture  area  ratio. 

The  effect  of  ramp  angle  on  additive  drag  is  presented  in  Figure  46.  The 
additive  drag  coefficient  increases  with  increasing  ramp  angle.  The  additive  drag  as 
defined  in  Equation  2-8  is  the  integral  of  the  static  pressure  over  the  projection  of  the 
streamtube  exterior.  As  demonstrated  in  Figure  45,  the  capture  area  decreases  with 
increasing  ramp  angle  or  the  axially  projected  area  of  the  streamtube  increases.  The 
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FIGURE  46.  Bffsct  of  Ramp  An|W  on  Intel  Additbt  Drag. 

increased  ramp  angle  causes  greater  deceleration  or  more  compression  of  the  How, 
increasing  the  static  pressure.  Thus,  both  increased  pressure  and  area  combine  for  an 
increase  in  additive  drag. 


Effect  of  Mach  Number 

These  same  results  have  been  replotted  to  demonstrate  the  effects  of  freestream 
Mach  number.  The  effect  of  Mach  number  on  inlet  capture  area  ratio  is  shown  lit 
Figure  47.  At  a constant  ramp  angle,  the  capture  area  ratio  increases  with  Much 
number.  As  the  Mach  number  increases,  the  shock  wave  angle  decreases,  and  the 
influence  of  the  ramp  is  exerted  over  a smaller  region.  At  Much  numbers  above  the 
design  or  “shock-on-lip”  Mash  number,  the  shock  wave  is  within  the  cowl,  and  there  is 
no  spillage  over  the  cowl,  only  sidespillage.  The  shapes  of  the  captured  streamtubes  are 
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FIGURE  47.  Effect  of  Mach  Number  on  Intel  Csptur*  Area  Ratio. 


presented  in  Figure  48  for  a ramp  angle  of  10  degrees.  Both  spillage  over  the  cowl 
i . and  sidespiliage  decrease  as  the  Mach  number  increases. 

I I 

| The  effect  of  Mach  number  on  additive  drag  is  shown  in  Figure  49.  Additive 

I drag  decreases  with  Mach  number.  The  increasing  capture  area  ratio  means  that  the 

axially  projected  area  of  the  streumtube  is  decreasing  with  Mach  number.  Although  the 
' j pressure  behind  the  shock  wave  increases  with  Mach  number,  the  pressure  coefficient 

i which  is  defined  as  the  difference  from  freestream  pressure  divided  by  the  freestream 

dynamic  pressure,  decreases  with  Mach  number.  Thus,  both  contributions  to  the 
additive  drag  coefficient  decrease  with  Mach  number  resulting  in  a decreasing  additive 
I drag  coefficient  (not  necessarily  decreaning  additive  drag). 


Effect  of  Sidewalls 


Sidewalls  are  often  added  to  two-dimensional  inlets  to  reduce  the  amount  of 
sidespillage.  An  analysis  was  performed  to  investigate  the  effects  of  various  sidewalls  on 
inlet  capture  area  ratio  and  additive  drag.  In  initialization  of  the  computer  program, 
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FIGURE  49.  Effort  of  Mach  Numbsr  on  Ink!  Addftto  Drag. 


points  arc  apeclfled  along  the  leading  edge  of  the  aideplate  for  tracing  back  to 
frecstream.  The  number  of  points,  spacing,  and  incremental  change  in  the  x direction 
are  determined  by  the  program  such  thut  new  points  are  introduced  at  each  data 
printing  interval,  and  the  planar  relation  of  each  printing  is  maintained.  As  in  the  case 
without  sidewalls,  two  points  are  specified  in  the  vicinity  of  the  side  edge  of  the 
ramp,  one  slightly  above  and  one  on  the  ramp  surface.  The  determination  of 
perturbation  velocities  and  the  streamline  trace  are  calculated  In  the  same  manner  os 
without  sidewalls.  The  analytical  model  requires  that  the  forward  portion  of  the  ramp 
is  not  bounded  by  the  sidewall.  A portion  of  the  ramp  is  necessary  to  establish  a 
conical  flow  regime  originating  from  the  comer  of  the  ramp  leading  edge.  The  model 
does  not  consider  the  comer  flow  resulting  from  the  intersection  of  the  sideplate  and 
the  ramp.  The  complicated  shock  intersections  and  corner  flow  occur  downstream  of 
the  cowl  and  sideplate  leading  edges  and  thus  have  no  effect  on  a streamline  trace 
upstream  to  determine  mau  capture  area  ratio  and  additive  drag.  A simple  sidewall 
with  straight  cutback,  as  illustrated  In  Figure  50,  was  assumed.  Other  sidewall  shapes 
could  be  accomodated  with  slight  changes  to  the  computer  program.  Various  ratios  of 
the  sideplate  length  to  iniet  length  were  assumed,  The  effect  of  sideplate  geometry  on 
inlet  capture  area  ratio  is  presented  in  Figure  51  for  the  case  of  a 10-degrec  ramp.  As 

87 


I 


NWC  T ’ 5951 


I I I 

OJB  0.00  0.70 

UWMTH  OP  mWMHUL  t^L) 

FIGURE  SI.  Effect  of  SWowalli  on  Intel  Captain  Aim  lUtto. 


■“‘•“"“inrin  aim 


88 


NWC  TP  595 1 

expected,  increasing  the  sidewall  area  reduces  the  spillage,  thus  increasing  the  capture 
area  ratio.  The  shapes  of  the  captured  streamtubes  are  illustrated  in  Figure  52  for 
Mach  number  of  2.5.  All  of  the  effect  of  sidewalls  is  shown  by  the  distortion  of  the 
side  of  the  streamtube.  Since  the  cowl  was  not  changed,  the  upper  portion,  which  was 
traced  from  the  cowl,  remains  undistorted.  As  the  sidewall  area  is  increased,  the 
sidespillage  is  reduced,  and  the  capture  area  increases. 

The  effect  of  aidewalla  on  additive  drag  ia  shown  in  Figure  53.  The  upper 
portion  of  the  streamtube  provides  moat  of  the  drag.  The  forces  on  the  side  changes 
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thin  value  only  slightly,  Both  poiitlve  and  negative  preiaure  coefficients  are  found  In 
the  aide  of  the  itreamtube.  Near  the  edge  of  the  ramp,  tho  lateral  perturbation! 
Increase,  resulting  In  local  pressures  less  than  freestream.  This  tends  to  decrease  the  net 
drag  of  the  itreamtube.  Depending  on  the  amount  of  lateral  spillage  and  the  area  of 
the  atreamtube,  the  contribution  of  sidewalk  may  1*  an  increase  or  decrease  in 
additive  drag.  This  is  indicated  in  Figure  S3.  At  a Mach  number  of  2.0,  the  presence 
of  sidewalls  reduces  the  lateral  spillage  and  increases  additive  drag.  At  the  higher  Mach 
numbers,  a minimum  drag  occurs.  In  all  cases,  the  effect  of  the  side  is  small  relative 
to  the  greater  contribution  of  the  upper  surface  of  the  itreamtube. 


Effect  of  Leading  Edge  Sweepback 

In  general,  two-dimensional  inlets  are  constructed  with  the  leading  edge 
perpendicular  to  the  aides,  l.e„  with  no  sweepback.  The  computer  program  was 


\ 
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employed  to  examine  the  effects  of  leading  edge  sweepback.  The  derivation  of  the 
perturbation  velocity  components  was  general,  not  limited  to  the  case  of  no 
sweepback,  An  inlet  with  sweepback  is  illustrated  in  Figure  54.  The  inlet  cowl  is  also 
swept  such  that  the  cowl  is  parallel  to  the  leading  edge.  To  consider  an  inlet  with 
sweepback,  two  sides  must  be  considered  separately  with  supplementary  sweepback 
angles.  The  individual  solutions  must  be  combined  to  obtain  the  solution  for  the  full 
inlet. 

The  effect  of  leading  edge  sweepback  on  the  capture  area  ratio  of  the  baseline 
10-degree  ramp  inlet  Is  illustrated  in  Figure  55,  There  is  no  change  in  the  capture  area 
ratio  at  any  of  the  Mach  numbers  considered  for  sweepback  angles  of  up  to  20 
degrees,  This  may  be  explained  by  examining  the  captured  streamtube  shown  in  Figure 
56.  This  was  determined  for  the  case  of  20-degree  sweepback  at  a Mach  number  of 
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2.5.  The  effect  of  sweepback  is  to  shift  the  flow  laterally.  Thus,  the  centerline  of  the 
• stream  tut  e is  displaced  from  the  centerline  of  the  inlet.  Comparing  the  itreamtubo 

with  sweepback  to  that  without  sweepback,  it  appears  that  lateral  displacements  on 
one  side  are  approximately  matched  by  similar  displacements  on  the  other  side.  Thus, 
little  effect  on  capture  area  ratio  is  anticipated.  Similarly,  sweepback  has  no  effeot  on 
additive  drag  as  shown  in  Figure  57.  Increases  in  drag  on  one  side  of  the  streamtube 
are  counteracted  by  decreases  on  the  other  side. 
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Chapter  6.  SUMMARY  AND  DISCUSSION  OF  APPLICABILITY 

The  problem  of  additive  drag  of  two-dimensional  inlets  has  been  examined.  The 
theory  of  homogeneous  conical  flows  as  described  by  Milne-Thomson  was  expanded  to 
determine  the  three-dimensional  perturbation  velocity  components  within  a region  of 
conical  flow  as  formed  by  the  edge  of  an  inlet.  A computer  program  was  written  to 
trace  streamlines  upstream  from  an  arbitrary  inlet  cowl  to  freestream  conditions.  This 
technique  is  used  to  determine  the  shape  and  area  of  the  captured  stream  tube.  The 
local  static  pressures,  as  computed  from  the  velocity  field,  are  numerically  integrated 
to  determine  the  inlet  additive  drag. 

The  computer  program  was  applied  to  a baseline  inlet  configuration  to  evaluate 
the  effects  of  various  inlet  design  parameters:  freestream  Mach  number,  ramp  angle, 
sidewalls,  and  leading  edge  sweepback.  Further  application  of  the  program  to  other 
configurations  such  as  both  sidewalls  and  sweepback  is  simple. 

Consideration  must  be  given  to  the  range  of  applicability  of  this  technique.  One 
limitation  is  the  small  perturbation  assumption  which  is  the  basis  of  the  method  of 
homogeneous  conical  flows.  It  is  assumed  that  the  perturbation  components  are  small 
relative  to  freestream  velocity  or  u,/V  < 1.  For  two-dimensional  flow, 

y ■ - o tan  u 
and 

V = « (6*1) 

For  Mach  numbers  greater  than  or  1.41,  the  tangent  of  the  Mach  angle  is 
less  than  unity  and  the  vertical  component  w is  the  limiting  condition. 

If,  for  example,  it  is  assumed  that  w/V  ■ 0.3  is  a reasonable  limit  of  the 
assumption  of  a small  perturbation  velocity,  the  theory  may  be  applied  to  ramp  angles 
of  up  to  17.2  degrees,  independent  of  Mach  number.  Other  choices  of  the  definition 
of  a “small  perturbation"  will  result  in  different  limiting  ramp  angles. 

An  alternate  approach  is  to  compare  the  perturbation  velocities  obtained  by 
oblique  shock  theory  and  by  linearized  theoiy  for  two-dimensional  flow.  These 
perturbation  velocities  in  the  vertical  and  axial  directions  are  compared  in  Figures  58 
and  59,  respectively.  As  has  been  previously  observed,  the  vertical  perturbation  velocity 
in  two-dimensional  linearized  theory  is  independent  of  Mach  number.  Since  the 
perturbation  velocities  by  both  theories  are  less  than  freestream  velocity,  the  magnitude 
of  the  differences  in  perturbation  velocity  divided  by  freestream  velocity  was 
considered  to  be  a valid  criterion  in  determining  the  range  of  applicability  of  the 
additive  drag  analysis.  The  range  is  presented  in  Figure  60  for  various  differences.  This 
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figure  was  controlled  by  the  longitudinal  perturbation  velocitiea.  The  vertical  velocitiea 
* demonstrated  smaller  differences  for  the  same  conditions.  As  can  be  seen  from  the 

| figure,  this  method  may  be  applied  for  reasonable  ranges  of  Mach  number  and  inlet 

| ramp  angle. 
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Appendix  A 

DERIVATION  OF  COMPATIBILITY  RELATIONSHIPS 


The  complex  velocity  components  U(e),  K(e),  and  W(e)  ore  not  independent 
since  their  real  parts,  u,  v,  and  w,  are  velocity  components  of  an  irrotational  flow. 
Since  u,  v,  and  w are  homogeneous  functions  of  degree  zero  in  coordinates  x,  y,  and 
z,  the  velocity  potential  from  which  they  are  derived  must  be  homogeneous  of  degree 
one,  Let  f(x,y,z)  be  a homogeneous  function  of  degree  n.  Euler’s  theorem  states  that 

xfx  + yfy  + zf2  - nf  (A-l) 

where  subscripts  indicate  partial  differentiation.  For  the  velocity  potential  4>,  n - 1 

x4>x  + y«hy  + zfcj  ■ <P  (A-2) 

Since 

9<t> 

‘h  * - xu  - yv  - zw  (A-3) 

Therefore 

d«h  * - udx  - vdy  - wdz  - xdu  - ydv  - zdw  (A-4) 

But  by  chain  rule  of  differentiation, 

d«l*  " 4»  dx  + 4>  dy  + 4>  dz  * - udx  - vdy  - wdz  (A-5) 

By  subtracting  from  Equation  A-4 

xdu  + ydv  + wdz  ■ 0 (A*6) 

Since  the  complex  components  u*,  v*,  and  w*  are  also  perturbation  components  of 
some  conical  flow,  it  follows  that 

xdu*  + ydv*  + wdz*  ■ 0 (A-7) 

By  multiplying  by  i and  adding  to  Equation  A-6 

x(du  + idu*)  + y(dv  + idv*)  + z(dw  + idw*)  ■ 0 (A-8) 


i 


NWC  TP  5951 


Since 

du  + idu*«d(/«  l/de 

where  ’ indicates  differentiation  with  respect  to  e,  Equation  A-8  becomes 

x£/*  + yV  + ziv'  -0 

Partial  differentiation  with  respect  to  x yields 
6/'+(xt/"  + yF"  d-zlP”)#,  -0 

or 

l//ex  --(xf/'  + y^'  + zr') 

Differentiating  Equation  A- 10  with  respect  to  y and  z yields 

mhUm  + yV"  + ZW") 
ex  ey  ez 

This  is  the  compatibility  equation  for  irrotational  flow.  It  may  be  reduced 
application.  From  the  definition  of  e, 

y + lz 

x tan  m ^ R2 

where 

R2  "(x2  tan2  fi-  y2  - z2)1/2 
By  logarithmic  differential  of  Equation  A-1 3 

de  dv+ldz  tan/idx  + dR2 
e " yTlz  x tan  m + R2 

Differentiating  Rj  ■ x2  tan2  p - y2  - z2  yield’ 

2R2dR2  ■ 2x  tan2  m dx  - 2ydy  - 2zdz 
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(A-1 1) 

(A-1 2) 
further  for 

(A- 13) 

(A- 14) 
(A-1 5) 
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By  substituting  into  Equation  A-t4 

de  dy  + idz  tan  M dx  + l/R2(x  tan2  m dx  - ydy  - zdz) 
e “ y +lz  " x tan  n + ft2 


T-Td»+'>y+'> 

By  collecting  terms  of  dx 


ex  tan  ju  dx  + l/Rj  (x  tan2  m dx) 
T11* " x tan  n + Rj 

ex  R2  tan  n + x tan2  /t 
T * " ~k2(x  tan  ft  + R2) 

tan  n (R2  + x tan  ft) 

" R2(x  tan  fi  +“ft2)_ 


e x taryt 

“e  “ "Rj 


or- 


e tan  #t  " R i 


By  collecting  terms  of  dy  in  Equation  A*  1 6 

ey  dy  >/R2(-ydy) 
e dy  ■ y + j7;  ~ x ton  m R"2~ 


Since 


Thus, 


y + iz 

•"xii lirirnq 

y + Iz  ■ «(x  tan  n + R2) 


| " e(x  tan  ft  + R25  + R2(x  tan  pi  + R2) 

I <y  R2  + y2 

, | « " «R2U  tan  n + R2!> 

( !; 

i!  ioi 


‘■VUt  . 


r ii  i n.  .iVtodUlU-i 


H JLiUn  l.».  ■ —*J  I . I . . - 


(A- 16) 
(A-l  7) 

(A- 18) 
(A-19) 


(A-20) 


(A*2I) 
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; 

(A-22)  : 
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or 


or 


R2  + y« 

ey  “ R2'(x  tiTu  + R"^ 


'y  i 

“R2?yT"“iq 

x tan  m + R2 

By  collecting  terms  of  dz  in  Equation  A- 16 


jdz- 


idz  _ gdz) 
y + iz  " x tan  n +1? j" 


Similar  to  the  y components, 


1 

IR2  + z<  Ttj 

x tan  m + 

Thus,  from  Equations  A-19,  A-23,  and  A-25 
«„  *y 

■ «'■  ■■.I.  i S fc..i  , ■ i 

- e tan  m R2  + y iR2  + z 
x tan*i  + R3  x tan #i  + k2 


From  Equation  A*  12 

t/MP  IV' 

«*  " ey  "•* 

*»  r2  +y* 

17  •„  * (x  tan  u + Rj)(-  « tan  *0 
Consider  each  term  separately 


R2 

i Term  1 " (x  tan  m + k2)(-  « tan  m) 

‘ 
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(A-23) 


(A-24) 


(A-25) 


(A-26) 
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(A-27) 
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From  Equation  B-5  in  Appendix  B 
R}  = x2  tan’  p ( 1 - r2 > 


Term  1 


- x tan  <i(l  - r2)l/2 


(x  tan  /i  + tan  /a  (1  - r2)l/:Xe  tan  /a) 


Term  1 


(1  - r2)>/2 


{1  +(1  - r2)^2)e  tan /a 


From  Equation  B-9 
1 - R2 

Thus, 


2 s/\  - r2 


I + \/l  - r2 


Term  l 


The  second  term  in  Equation  A*27  is 

Term  2 ■ ^ tun  ju  + &2X-  e tan  y) 
From  Equution  B-2 


a m r 0 ■ 

" r v x tan  /a 


y ■(rcosOXxtanja) 

By  substituting  for  y and  for  R2  from  Equution  B-5 


Term  2 


g(r  cos  OXx  tan  n) 


(x  tania  + x tan  m(  1 - r2)l/2)[-  e tan /a) 


Term  2 


- ercosfl 


e tan  ju  (1  +(l  - r2)1^2) 
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From  Equation  B-8 


r. — 2 i - r2 

1 - r ■ r 

I + R2 


and  from  Equation  B-IO 


1 + R2 

Equation  A*29  becomes 


- e cos  9 


Term  2 


Term  2 


e cos  9 2R 


• tanu  ,+JLlRI  eUnM(l*R2  + l-R2) 
1 + R* 


2cR  cos  $ 


2e  tan  fi 


Combining  the  terms  A-28  and  A-30,  Equation  A-27  becomes 
V _gy  _ - (t  - R2)  2*R cos 6 


? - *y  _ R2  - 2«R  cos  9 - 1 


2«  tan  n 


e - Re1®  ■ R cos  0 + Ri  sin  0 
e/R  ■ cos  0 + 1 sin  9 
R/e  * e“ 19  ■ cos  9 - i sin  9 
e/R  + R/e  * 2 cos  8 
e2  + R2  ■ 2Re  cos  8 
e2  ■ 2Re  cos  9 - R2 
-0  +e2)«-  1 - 2Re  cos  0 + R2) 

By  substituting  into  Equation  A-31 


) 

1 

It! 
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or 

W a v* 

2e  tan  n (1  +e2) 

In  a manner  similar  to  the  above  deviation  for  V*,  from  Equations  A-12  and 
(f  W' 

€x"ei 

et  iRj  + ze 

Z7-"  e^m  (x  tan  n + R2)(-  e tan  n) 


as  before,  In  Equation  A-28,  the  first  term  becomes 

. - id  - R2) 

Term  1 - 

Since  s*  ■ r sin  0 * z/x  tan  n,  the  second  term  of  Equation  A>34  becomes 


u 


Term  2 ■ <x  tan /i  + Rj)(-  e tan  ,u) 


T,m  , _ - _2er  sin  0 
Term  2 tan /T 

and  Equation  A-39  becomes 


W'  h - i(  1 - R2)  - 2er  sin  0 


{7  e 


2c  tan  /i 


Again, 

e/R  ■ cos  0 + I sin  6 
R/e  ■ e"1®  ■ cos  6 - i sin  0 
e/R  - R/e  - 2i  sin  6 

e2  «2Rei  sin  0 + R2 
1 - e2  - 1 - 2Rel  sin  9 - R2 
- id  - e2)--  1(1  - R2)-  2Re  sin  9 
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A-26 

(A-34) 
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From  Equation  A-34 

W et  - i(l  - e2) 

TFm  ex  “ ' 2e  tan  n 

Combining  with  Equation  A-33 

V V*  m W' 

it  tan  m _ (i  + e2)  " _ i(j  _ e2) 

which  are  the  compatibility  relationships  for  irrotational  flow. 


(A-36) 
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Appendix  B 

DERIVATION  OF  SIMPLIFYING  RELATIONSHIPS 
FOR  THE  BUSEMANN  TRANSFORMATION 

The  following  relationships  are  useful  in  deriving  the  compatibility  relationships 
described  in  Appendix  A,  Consider  Busemann’s  complex  variable 


_ y + iz 
6 ~ x tan  p + ft2 


(B-l) 


where  R2  = (x2  tan2  p - y2  - z2)1  ^ 

The  x,  y,  and  z symbols  are  Cartesian  coordinates  and  p is  the  Mach  angle 
sin"1  1/M.  The  dimensionless  variables  r,  s,  and  s*  may  be  employed  in  mapping  from 
the  physical  plane  to  the  e plane,  where 


s ■ r cos  6 = 


x tan  p 


s*  “ r sin  9 


x tan  pi 
s + is* 


1 + (1  - s2  - s*2)1^2 


, . ■<<*» " + . Rei» 

1 +(1  - r2)1^2 


(B-2) 


(B-3) 


where 


R. I 1 -U-  V-t'l 

1 +(l-r »)'/»  ' 

Proof  1 : Rjj  ■ x2  tan2  p ( 1 - r2) 

R^  - x2  tan2  p - y2  - z2  - x2  tan2  p - s2  tan2  p - s*2  tan2  p 

■ x2  tar.2  p (1  - s2  - s*2)  “ x2  tan2  p (1  - r2  cos2  0 - r2  sin2  9 ) 
Rj  ■ x2  tan2  p (1  - r2) 


(B-4) 


(B-5) 
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Proof  2:  >/(!-  ^"TTp- 


from  Equation  B*4 


l-O-r2)!/2 


1 - 2(1-  r2)1/2  + 1 - r 


.2  2-2(1-  r2)1/2-  r2 


1 - R2  ■ 1 - 


2-  2(l-r 2)112- 


! . m£ll±Mz £)-!!  * j.  , 2f2  " 2 + 2(1  - r2)1 1 2 
r2  r 


h-r2,  £±lrJlL£){/2.:,S. 

r 


1-R2  B 2r*-  2 + 2(1  - r2)1^2  a r2-  1 + (1-  rV 
1+R2  ”""2- 2(i  -r2),/2/r2  1-d-r2)1/2 

^ [r2  - 1 + (1  - i2)1^2!!!  + (1  - r2)1^2] 

1 - (1  - r2) 


l+R2  v 


"■•WWMlMSJj) 
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Proof  3:  2 - 1 - R2 

1 +n/1-  r2 


From  Equation  B*6 


, _ p2  T + 2(1  “ t2)11’  _ [2r*  - 2 + 2(1  - r2!1^2]  [1 + (1  - r2)1/21 
r2  r^l+U-r2)1'2] 

_ gr2  ~ 2 2(1  ~ r2)1^2  '+  2r2fl  - r2!1'2  - 2(1  - r2)1'2  + 2(1  - r) 

r2 ! 1 +(1  - r2)'/2) 

— aAl-^/a 

r2 [1  + (1-  r2)1/2] 


« o'1  2\/l  - r2 

1 - R2  y~T*m,  iii 


1 4 « r2 


Proof  4:  r* 


1 + R2 


From  Equation  B-4 


1 ■ 1 1 - r2)1^2 


From  Equation  B-7 


r2  I3 


By  lubitifution  for  R 


l + R2-f 


1 + R2 


(B-10) 
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Appendix  C 

LISTING  AND  SAMPLE  PRINTOUT  OF  COMPUTER  PROGRAM 


9 


Program  Ustini 

C THI*  PROORAM  TRACES  THE  STREAMLINES  FROM  THE  COWL 

C Of  AN  INLET  USINO  THE  HETHOO  Of  HOHOOCNEOUS  CONICAL 
C FLOHS  - MILNE -THOMSON 

COMPLEX  CMS. ANU.AI ,SI .CU.CV.CHI .CHF.C FNU 
COMPLEX  CVELI .CVELF.CVEL3 
COMPLEX  Tl .TF.T3.TH.TS.TS 

01 MENS I ON  XI 300.301 .Yl 300, SO) ,31300.80) ,PI 300, SO) 
NRITEIS, ISO) 

ISO  FORMAT  I I HI  I 
C INITIALISE  PROSLEM 

CRAP  - 3. IHISSISSH/ttO. 

OS  REAOIO. I 10)  AMACH 
I 1 0 FORMAT  I OF  1 0 . 3 ) 

IFIAMACH.LT. 0.0)  00  TO  OS 
HR  I TE I S , II  I I AMACH 

III  FORMAT!/, 06H  FREEST REAM  MACH  NUMBER  • ,FI0.3> 

VAST AX-SORT  I I .0* AMACH* AMACH/ I I . ♦. 1 'AMACH* AMACH I ) 
PIERO-  I ./HI  . ♦ . 3*  AMACH*  AMACH)  * *3 , S ) 

OICRO-. 7 •AMACH* AHACH’PZERO 
AMU-ASINI I . /AMACH) 

TANMU-TANIAMU) 

UINF-I .0 
VINF-0 . 0 

C OttCRIsr  INLET  ANO  INITIALISE  ORIO 

READ  I B , I 10 1 ALPHA, XCOHL , ZCOML , XSIOE , HIOTM.DCLTA 

MIOTH-MIOTM/O.O 

SB IDC-XS IOC -TANI ALPMA-CRAD) 

ARCF -H I OTH*  ZCOHL 

NSl OE-B  + 10,0*1 XCOHL -XS I OE I /XCOHL 
NL INC- 1 1 ♦NSIOE 
01- 1 ICOHL-SS IDE  I /NS IOC 
ISO  FORMAT IfllO.FIO.H) 

URITCIS, I IF  I ALPHA, XCOHL, ICOHL , XSIOE, SSIOC 
lit  FORMAT (•  INLCT  RAMP  ANOLC  - '.FIO.S./, 

I'  COHL  COORDS  > •.■FiO. 3./,'  SIDE  COORDS  - * .0FI0.3) 
URITCIS, 1 1 3 IHI OTH, DELTA 
113  FORMAT  I ISH  INLCT  HIOTH  > .FIO.S,/. 

IFIH  INLCT  SHEEP  ANOLC  - , FIO.S) 

URITCIS, I IS) 

MlNF--ALPHA<CRAO 

IF IOCLT A , EO . 90 . 0)  00  TO  SO 

C3- 1 . 0/  I TAIIIOELTA'CRAO)  I 

COSLAM- -TANMU/ I TANI DELTA *CRAO> I 

AL  MOO  A -AC  OS  I COSLAM  I 

TMETA-.S-ALMOOA 

00  TO  SI 

50  THE TA-HB.O -CHAO 
ALMODA-SO . 0 *CRAO 
C3-0.0 

51  OI-COSITHCTA) 


\ 


III 


TOG1DZN0  *JOI  KUJC-MOT 
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DB-t INI THCTA) 

Al-CMTLXIOI ,DBI  r 

•l-CONJOIAI  I 

MICRO— ALOMA>CNAO>TANHU/B INI ALHBOAI 
OVaMIOTH/IO. 

1*1 

MRITCIO.ISOI 

IS*  rORHATlNX.  I Ml  ,NX.  IHJ.OX.OHXI  I .Jl  .HX.OHVI  I ,JI  .*»*, 

I SHU  I I ,J1I 
00  BO  J-l.tl 
XI I . JI-XOOML 
VI I ,Jla|J-lll »0V 
XI  I , JlalCOHi. 

■o  mritcio.iinii  .j.xi i , ji  ,vi  i ,ji  .xn  ,ji 

ll<*  RORNAT  IBIB.BX.iriO.ll 

iriXCOWL.NC .XOIDCIOO  TO  *1 
C NO  BIOCRLATC 
MR  I TC IB ■ 111 t 

IBI  rORMAT I / i I MM  NO  BIOCRLATC  > 

KBIOC«0 
NXX-BC I 
JX-0 

OX-XOOML/BOO. 

NO- ML  INC -l 
00  OB  J-IB.NO 
XI  I , Jl aXCOML 
VI  I , Jl-0.0 

z 1 1 ■ j i azeoML - 1 j- i it  *ox 

•»  MR  I TC 10. 1 INI I , J.XI I , Jl ,VII ,JI .XI I .Jl 
IX-I 

00  VO  IN 

C INLCT  WITH  BIOCRLATC 
tl  OXalXCOHL-XBIOC WNBIOC 
HRITCIO, ibbi 

IBB  roRMATI/, INN  BIOCRLATC  > 

KBIOCal 
JX*  I 

NXX-IO*XCOML/OX  • I 
NO  ■ NL  INC*  I 
00  C)  Jail, NO 
l-IO*IJ-l !>♦! 

xii  ,ji«xcoml-i j-i  n*ox 

VI  I , JI-0.0 

XI  I t Jl aZCOML* I J- III *DC 
OS  MRITCIB.I INI I iJiXII , J* , VI I ,JI .XII  . Jl 
I X—  I ♦ 1 0 
NO- II 
ox-ox/ to. 

BN  JaNLIHC 
l-IX 

XI I . JI-XBI0C 
Vll.JiaO.O 
III . JlaZBIDC* . BO*DZ 

MRITCIO. I INI  I ,J,MI I . Jl .VI I , Jl ,11 1 , Jl 

J*M.  INC*  I * 

XI  I , Jl «XBIOC 
VI  I .JI*-.BO*OX 

►> 
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•1 

f 


I; 

l' 


ir 


» 


\ 


e 


e 


c 


e 


c 


c 


c 


c 


z« i , ji*zaioc 


hritcis.iioi 

hritcis.iisi 

C VALUATION  OP  CONSTANTS 

C I "COSI alpha»crao> 

Ct»SIN< ALPHA • CP AO  I 

CU« -UZCPO 'CHPLX 1 0 • 0 > I 01/1.  INISMSSN 
CV«-CU/ I X*  T ANHUI 

CHI a*UlKRO*CHPLXI 1 . 0 . 0. 0 1 / 1 1 . 0»l. INISBtf JN»T ANHUI 

CHXa-S . 0 'UZCRO'O I *DS STANHU 

STAC AML  I NC  TRACK 

NO’NL INC A I - JX*  I NS  I DC A 1 1 

KX>NO*IO 

00  II  |a| iNXX 

00  30  J«l,NO 
R0TA7C  COORDINATCS 
XR-CI *X< l,JI»CI'III , Jl 
zr»-c«»xi i ,ji*ci>z< i ,j» 

YRa  Y I 1 ,J) 

If IJ.CO.NLINC*! IZR-C.O 
IPIXP.LC.O.OI  00  TO  IN 
CMCCKalXR»XR»TANHU*TANHU-VR»VR-XR»IRI 
ir  ichcck.lt . o • oi oo  to  11 

POINT  IS  HITHIN  NACH  CONK.  USC  CONICAL  PLOH 

Rt-SORTICHCCKI 

DCNOH  •XR*TANHU*R« 

CPS*CHPLXI YP/OCNOM.ZP/DCNOHI 
ANUaCSQRT ICPSI 

CONICAL  PCRTORRAT ION  COHPONCNTS 
KTVPC«I 

CPNUal I Al 4ANU) • IRI »ANUI 1/ 1 1 A I -ANUI • IS I -ANU I I 
CVCLl«CU*CLOOICPNU» 

UXaRCALICVCLI » 

' T|a«./IAI*ANU!*t.MBI*ANUI 
T»a«,  >AI  •ANU*I.  ••  I • ANU 

Tla| I . /( AI*AI I I •CLOOl I AI»ANUI/ANUI»< I . / <01 *01  I * 

1 *CLOO  HOI  a ANUI  / ANUI 

TN-I I . / I At • A I > I *CL00l IAI-ANU) /ANUI* 1 I . / ISI *01 1 1 
I *CLOO I 101 "ANUI /ANUI 

TS.AI  • A I *CLOO(  A I a.ANUI  I *CLOOIS  I ♦ ANUI 

TS*AI *AI *CLOOI A I -ANU I »SI *CLOOISI - ANUI 
CVCL«»CV»I-TI-T»*T1“TN*TS-TSI 
UV«RC AL I CVCLO * 

CVCLI«CHI»I-TI*TR*TJ-TN-TSyTSI*CM# 

UXaRCALICVCLI! 

00  TO  IB 

THO  OiHKNSIONAL  PLOH 
IR  TANA  ■ ZR/XR 

IP  I TANA.OK . TAMHUIOO  TO  IN 

THO  PIHCNIIONAL  I PLANAR  I PLOH 

KTVPCaf 

UXaUXCRO 

UVaVXKROaCI 

UZ*  ALPHA  *CRAO 

OO  TO  M 

UNO  I STUROCO  PLOH 
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t.  ’ 


IS  KTVPCa  I 

ux-o.a 
uv«o.o 
uz«o . 0 

C STM  AML  INC  TNACC  TO  UPOTNCAH  POSITION 

M VCLXaUINPaUX 
VCLY«VINP*UV 
VKLfaHINPaUI 
OV«OX*VCLV/VCLX 
OI«OX>VCLZ/VCLX 

VVaHOftTlVCLX*VKLX«VCLV*VCLV*VKLI*VKLI» 

VP»VV<>VA«TAN 

PP-t I . -VP*VP/0.  |M|,| 

PI  I , 4KPP-PZCN0 

50  XTaXN-OX 
VTavo-ov 
it-in-oi 

XNaXT 
VN-VT 
II  IP-IT 

lie  POHMATIOIO.HP  10. at 
XI l«l .JI-CI*XT-GO*IT 
II  l«l  ,4»aC0>XT*CI*ZT 
VI 1*1 ,41-VT 
ao  CONTINUC 

I IB  PONNATiaiB.7PIO.HI 
KX>HX*JX 
NO-KX/IO 

IPINO.NC .NL  INC  100  TO  SI 

NO-NOPJX 

KX-NX»IO 

jx-a 

51  CONTINUE 
111  PONNATIBIBI 

■ JX-K»IOC 

NO-NLINC*  I -JX»  (N*IO«»l  » 

00  71  I • I ,NKX,  10 
HNITKIO.ISOI 

ISO  PONMATIHX, INI ,HX , IHJ , HX ,OHX I I .41 .SM.OHVI I ,41 ,HX. 

I ONI  I I ,41  .HX.OHPI  1.4)1 
00  70  J-I.ND 

70  HNITCIO.IBOII  , J.XI  I ,41  .VI I .01  ,11 1 ,41  .Nil. 41 
N0-N0*4X 

IP INO.NK.NL INC IOO  TO  71 

ND«NO*JM 

<M*0 

71  MNITKIO. 1 101 
1 10  PONNATt | HOI 

C OTMANTUOC  AN* A CALCULATION 
ANCA  * 0.0 
00  HO  4" I .NL INC 

DANE A» . 0* I VINNX.4* I I - VINXH.4) I *ICINKK.4v 1 1 *1 INXX .41 1 
HO  ANCA  a AM  A *0  AM  A 
AOACaANCA/ANKP 
HNITClO.  I SOI  AM  A.  AO  AC 

ISO  PONHAT I / , I OH  OTXCANTUOC  ANCA  « .PIO.S./.ON  AO/AO  a 
I .PIO.OI 
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ADO) T I VC  ORAO  CALCULATION 

DRAO-O . 0 

ATOT-O.O 

DTOT-O.O 

NO-NXX-I 

JXK-  1 1 

00  *1  l-l .NO 
00  00  J-l , 10 

MAN- 1 PI  I ,Jt*PI 1*1 ,JI*PI I ,J*I l*Pl 1*1 ,J*I I l/N. 
A0l-.0*IVtl,J*ll-VII.JII*lHI.J*|l*Hl«JII 
ADI-.t*  I VI  l*l,J*ll“VII,J*tll*llll'»l«J*|t*HI,J*lll 
A01-.0* I VI l«l .JI-VI I ♦ I ,4-1 M *IIM*I . JI*II 1*1 .0*1 n 
AOH - i 0 * I V I I ,J»-VI 1*1 ,JI 1*111 I ,JI*II 1*1 ,Jt I 
AONAO-AOI *A0I*A01* AOR 
DTOT-DTOT*ADRAO 
ATOT-ATOT*AORAO 
DORAO-POAR*AORAO 
•0  DRAO-ORAO*OORAO 
■I  OTOT-O.O 

MR I TC 10. 1ST I 

1ST  FORMAT  I/, I RM  URRKR  OUNFAOKI 
CORAO-DRAO/ I ARK F •OZCROI 
HRITKIO, III  IORAO.CORAO.ATOT 
00  01  l-l .NO. 10 
00  00  J-l I .NLINC 
IF  I KO I DC . CO . 0 >00  TO 

IFI I . OC.NLINC-IOIOO  TO  00 

IFIJ. LT .JUKI  00  TO  00 
IFIJ.OT.JMKI  00  TO  07 

ROAR- I R t |,J>*Pil*lO,J)*PII*IO.J*l>l/a,0 
AO  I -.0*1  VI  1*10,  J*l  I -VI  I , JM-III  1*10, J* I 1*11  I . Jl  I 
ADI- > 0* I V I 1*10. JI-VI 1*10. J-l I 1*1X1 l*IO.JI*II 1*10, J* III 
A01-  . 0*  I V I I, Jl-V  11*10.  JIHIXII.JI  *111*10, Jll 

aorao-aoi*aoi*aoi 

00  TO  00 

00  ROAR- 1 R I | ,JI«RI | • 10 ■ Jl *RI I .J*l I *Rl 1*10. J* I I I FR . 0 
ADI -.0*1 VI  I , J*l I -VI I , Jl t *1*1  I . J*l 1**1  I , Jll 
ADI-.t* IVI 1*10, J* I • -VI  I . J*l I I • HI  I* 10. J* 1 1*11 1 , J*l I • 
AO!” . 0* I V I 1*10, JI-VI 1*10, J*l 1 1*1 II 1*10, Jl*«l 1*10, J*l 1 1 
APR-.**  t VI  I. JI-VI|*I0.JI 1*1111, Jl *111*10. Jll 

aprao-aoi*aoi*adi*aor 

SO  OTOT-OTOT * ADR AO 
ATOT-ATOTAAMRAO 
OORAO-ROAR* AORAO 

00  DRAO-ORAO*ODRAO 
07  QONTINUC 

JRH-JKM*I 

01  OTOT-O.O 
HRITKIO.ISOI 

■110  FORMAT  1 1 , 1 7H  TOTAL  OTRKAMTUM I 
COR AO-OR AO/ I ARIF 'OIKROI 
HRITKIO. Ill IDRAO.CORAO.ATOT 
111  FORMATION  OR AO  - .FIO.O./.OM  CORAO  - .FIO.O./, 

HIM  TURK  OURFACK  ARCA  • ,F 10.01 
MRITCIO, 1 101 
00  TO  00 
00  HRITKIO, HOI 
STOP 
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Sample  Printout 


rnCCOTRCAH  HACH  MU HO 14  - 

0.000 

1 INLET  NAHA  ANOLE  ■ 

10.000 

COM.  COMOO  ■ 

10.000 

0.000 

} 0IOC  COMDO  ■ 

10.000 

1 .701 

INLET  WIDTH  * 

0.000 

! INLET  IMCKP  ANOLC 

i ■ 

. 00 

000 

i 

5 1 J 

XI 1 .Jl 

VI  1 .Jl 

*11. Jl 

1 I 1 

to.ooo 

-0.000 

0.000 

i 1 * 

10.000 

•4.000 

0.000 

j : ii 

10.000 

•4.000 

0.000 

j 1 14 

10.000 

-1.000 

0.000 

1 0 

10.000 

-1.000 

0.000 

10 

10.000 

-0.000 

0.000 

i ; 17 

10.000 

-0.000 

0.000 

1 0 

10.000 

-1 .000 

0.000 

1 0 

10.000 

-1.000 

0.000 

1 10 

10.000 

-.000 

0.000 

1 II 

10.000 

.000 

0.000 

NO  OIOCALATE 

1 10 

10.000 

.000 

4.101 

1 II 

10.000 

.000 

1.700 

1 14 

10.000 

.000 

1.000 

1 IO 

10.000 

.000 

0.411 

1 10 

10.000 

.000 

1 .001 

1 17 

10.000 

-.1*0 

1 .701 

XI  i , Jl 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 

10.00000 


VI  I .Jl 

•0.00000 
•V. 00000 
-4.00000 
-1.00000 
-1.00000 
-0.00000 
-a. ooooo 

•I .00000 
•I .00000 
-.00000 
.00000 
.00000 
. ooooo 
.00000 
.00000 
.00000 
-.  tor** 


XI  I .Jl 
0.00000 
0.00000 
o. ooooo 
0.00000 
Q. 00000 

o. ooooo 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
4 . 10000 
1.70011 
1.00700 
0.41000 
I .00074 
1 . 70107 


J XII, Jl  VII, Jl  *11. jl 
I 0.00700  -0.00000  4.01110 
0 0.00700  -4.00000  4.01110 


01 1 .jl 

. 04000 
. 04000 
. 04000 
. 04000 
.01001 
. ooooo 

. 00401 
.00040 
.01000 
.01100 
.01011 
.00700 
.00100 
-.00011 
-.01400 
- .00414 
-.00100 


Oil  ,JI 

.04000 

. 0400C 
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*.*£7*0 

-H. 00000 

H. 91119 

H 

9.00700 

-3.00000 

H.9I3I* 

0 

9 . BOBO* 

-3.011*9 

H.MH39 

• 

O.BOHOO 

-*.91*79 

H.9330* 

7 

• . BOOM 

-*.o*ihh 

H . 9H 1 39 

• 

• . BO  1 *7 

-I.0MH9 

H . *H*0* 

• 

9.HUI3 

- 1 . 0*7*H 

H . ••■03 

10 

*.ho«** 

-.*•7*9 

H . MH07 

1 1 

b.hotho 

- . 0**37 

H . 970** 

I* 

a.HMBI 

- > 03BIH 

H . ]*■•■ 

13 

I.HH3I 

- . OHBIO 

3.IMIB 

IH 

*.H*3H0 

-. 09*33 

3 . OB  1 *B 

IB 

*.H**HB 

- . 00013 

* . H*007 

I* 

*.H*»I* 

-.*0017 

1 . BH997 

17 

9.007** 

-.3*»aH 

1 . B70HB 

J 

KII.JI 

VI  1 ,Jt 

211 .0) 

1 

*.01*1* 

-*.00000 

H.9M3B 

• 

9.01*19 

-H .00000 

H . **03* 

3 

*.OI*l* 

-H.OOOOO 

H . 9**3B 

H 

*.01*1* 

-3.90000 

H . 9*030 

• 

*.01*93 

-s.oia*3 

H.B19I9 

* 

*.00*93 

-IBID* 

H . 9MH  1 

7 

*.00*9* 

-*.0H*7* 

H . *793H 

• 

*.0030* 

-I.BH9M 

H . 9*9 1 9 

• 

*.00031 

- 1 . 0**B* 

H.9I07* 

10 

• . 9979H 

-e*a*s 

H . **B*7 

1 1 

*. 99*10 

- * 0H**« 

H . 91*00 

1* 

• . **3HS 

- . 097H9 

H . 30*3* 

IS 

•••117 

-.0**** 

3. 9*7*7 

IH 

*.*•77* 

-. 1*10* 

3.0H009 

10 

*.9013* 

- . I*H00 

*.H**7I 

I* 

• . 970** 

- . 30HH9 

1 . 9197* 

17 

*.OI*l* 

-.0*909 

1 . 0*90* 

<J 

Ml  1 ,Jt 

VI  1 ,JI 

III ,JI 

1 

■ . B**79 

-9.00000 

H . 73991 

• 

•01*79 

-H.BOOOO 

H. 7 39*3 

1 

* . **«79 

-H.OOOOO 

H . 73993 

H 

*.•••79 

-3.00000 

H . 71003 

• 

*.*ao** 

-3.  DIMS 

H . 7*137 

• 

• • • 1 *HB 

•*.* 39*9 

H. 701  13 

7 

*.*0997 

-*.09707 

H .*1**0 

• 

* . *0*90 

-1  . B*7*C 

H . «37*H 

• 

a.Boias 

-1  .07339 

H.90I7* 

10 

*.H*7I* 

-.073*0 

H . *9*07 

1 1 

*.H*3*9 

• . 0700* 

H * *0*7* 

1* 

*.H*07* 

- . 0***7 

H .*73*1 

13 

• . H9709 

-.1*911 

3.9HH0H 

IH 

• . H***7 

-.17919 

3.0*39* 

I* 

*.H7BI* 

-.•9*39 

• . H*0H7 

10 

■ . H73M 

- . H9339 

1 .*000* 

17 
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171 

10 

1 .*17*0 

-.00810 

* . 78**0 

-.007** 

171 

1 1 

1 .*0*07 

-.10883 

*.00037 

-.007** 

171 

18 

1 .*10*8 

-.101** 

* . 00888 

- . 007** 

171 

13 

1 .*8881 

-.88081 

3.3808* 

-.007** 

• 71 

1* 

1 .*381 1 

-.*88*8 

1.88888 

- . 007** 

171 

10 

1 .*0188 

- . 78878 

1 .8*33* 

-.007** 

171 

IS 

1 .08818 

-1 .18370 

.80880 

. 0*808 

171 

17 

1 .08818 

-1.30180 

.88788 

.0*888 

1 

J 

Ml  1 . Jl 

VI  1 , Jl 

11 1 .Jl 

81 1 .Jl 

101 

1 

.80808 

•0 . 00008 

*.38110 

-.007** 

181 

8 

.80808 

-* . 08000 

*.38110 

-.007** 

181 

a 

.88808 

-*. COOPS 

*.38110 

-.007** 

181 

* 

.08808 

-3.00008 

*.38110 

-.087** 
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181 

8 

.88777 

-3.01883 

ISI 

8 

.88818 

-8 . 83088 

181 

7 

.07818 

-8.08877 

181 

8 

. 83887 

-1 .88873 

181 

8 

.88881 

- 1 . 1 0088 

181 

10 

.80880 

-.80818 

ISI 

II 

.88881 

-.10083 

ISI 

18 

■80783 

18177 

ISI 

13 

.81838 

-.88881 

181 

17 

.88788 

- . 78878 

ISI 

18 

.87370 

- . 78*78 

181 

18 

1 .08388 

-1 . 18370 

ISI 

17 

1 . 13878 

•1 .38180 

1 

J 

XII. *1 

VII.JI 

181 

1 

.78088 

-8.00000 

181 

8 

.78088 

-7.80000 

181 

3 

.78080 

-7,00000 

181 

7 

.78088 

-3.80000 

181 

8 

.78718 

•3.01883 

ISI 

8 

.78080 

-8.83888 

181 

7 

.77088 

-9 . 08877 

ISI 

8 

.78808 

-1.88873 

ISI 

8 

.71808 

-1 . 10088 

181 

10 

.70837 

-.80818 

ISI 

1 1 

.38888 

-. 10083 

ISI 

18 

.70037 

-.18177 

ISI 

IS 

.70778 

-.88881 

181 

17 

.71700 

- . 78878 

181 

18 

.73817 

- . 78*70 

181 

18 

.81808 

-1 . 18370 

101 

17 

.87731 

-1.38180 

1 

J 

XI 1 vJ> 

VI 1 .Jl 

801 

1 

• 

o 

* 

8 

o 

-8.00000 

801 

8 

- ■ 07880 

-7.80000 

801 

3 

-.07880 

-7.00000 

801 

7 

- . 07880 

-3.80000 

801 

8 

-.00038 

-3.01883 

801 

* 

-.08888 

-8.83808 

801 

7 

-.08887 

-8.08877 

801 

8 

-.07877 

- 1 .88873 

801 

8 

- . 08880 

- 1 . 1 0088 

801 

10 

-.10888 

-.80818 

881 

II 

-.11881 

• . 1 0083 

801 

18 

-.10718 

-.  18177 

SOI 

13 

-.08877 

-.88881 

801 

17 

-.08088 

-.78878 

801 

18 

-.07171 

- . 78878 

801 

18 

.00883 

-1  . 18370 

801 

17 

. 18181 

-1.38180 

7.71*83 
H . 77880 
7 . SOtSt 
7.87108 
7.08178 
7 • 71MB 
7.78070 

H . 0B7S7 
3.30788 
B.BBS1B 

I . 07870 
.B7I IB 
. BOOHS 


31  I .J> 

7.38081 

7. 3808 1 

7.38081 

7.38081 
7.711  as 
7.77800 
7.80073 
**.87100 
*4,88087 
7 . 7BBBB 
**  . 7BBS0 
*4  . 08878 
3 . 38707 
B.  8878 1 
I . 8*4  IBB 
.87030 
. I 1381 


XI  I ,J» 

*4 . 3883  3 
7 . 38833 
>«  . 38833 
*1 . 38833 
7,71  078 
7.7701  I 
*4 .80*488 
*4.87018 
*4 , 8*4888 
7,7*178 
7 . 78778 
7 . 0888 1 
1.38818 
8.88888 
1 .87088 
.88871 
. 08878 


‘.00777 
. 00777 
. 09777 
. 00777 
. 00777 
. 00777 
. 00777 
. 00777 
. 00777 
.00777 
.00777 
. 00777 
. 07888 


81  I , Jl 
.00777 
.00777 
.00777 
. 00777 
.00777 
. 00777 
. 00777 
.00777 
. 00777 
. 00777 
. 00777 
. 00777 
.00777 
. 00777 
. 00777 
.00777 
. 07088 


81  I . Jl 
- . 00777 
-.00777 
-.00777 
-.00777 
- . 00777 
- . 00777 
-.00777 
- . 00777 
- . 00777 
- . 00777 
• . 00777 
• . 00777 
-.00777 
-.00777 
-.00777 
- . 00777 
. 07888 
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Appendix  D 

DERIVATIONS  OF  TWO-DIMENSIONAL  PERTURBATION  VELOCITIES 


This  appendix  preienti  the  derivation  of  the  two-dimensional  perturbation 
velocities  to  be  used  as  boundary  conditions  for  conical  flow  analysis.  First,  the 
perturbation  velocities  are  derived  for  on  airfoil  without  sweepback.  These  results  are 
employed  to  determine  the  perturbation  velocity  components  for  an  airfoil  with 
sweepback, 


AIRFOIL  WITHOUT  SWEEPBACK 

The  small  perturbation  equation  of  velocity  potential  ♦ for  supersonic  flow  in 
two  dimensions  is: 


92<t>  2 

This  has  the  general  solution 

<t  - f,(z  - x tan  Mq)  + f2(*  + x tan  p0) 
where  f,  and  f2  are  arbitrary  functions. 

Now  consider  a velocity  potential 

4» « V*  4 Va(z  + x tan  m0) 


(D-l) 


(D-2) 


For  ee  - 0,  u ■ - 94»/9x  - - V which  is  a uniform  velocity  V in  the  negative  x 
direction.  Note  that  in  Milne-Thomson’s  terminology,  u(  ■ - 94»/9xt  which  is  opposite 
in  sign  to  most  terminologies. 

For  a small  but  non-zero  a,  the  velocity  components  ate  -V(l  + atanu),  -Va. 
Physically,  this  is  equivalent  to  flow  turned  through  an  angle  a.  The  uniform  velocity 
V is  separated  from  the  perturbed  flow  by  a characteristic  line  inclined  at  the  Mach 
angle  ^Q.  Thus,  the  perturbation  components  may  be  written  as 
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AIRFOIL  WITH  SWEEPBACK 


Now,  consider  a flat  plate  of  infinite  span  with  sweepback  as  shown  in  Figure 
D-l.  The  leading  edge  ia  swept  back  with  an  angle  8 between  the  leading  edge  and  an 
axis  parallel  to  the  freestream  velocity  vector  V.  The  airfoil  is  at  an  angle  of  attack,  a, 


At  point  O on  the  leading  edge,  define  two  planes,  OAB  in  the  plane  of  the 
airfoil  and  OCD  the  freestream  velocity  vector  V and  the  leading  edge.  From  this, 
draw  three  more  planes;  plane  OAD  through  V and  perpendicular  to  the  airfoil,  plane 
OCB  perpendicular  to  the  leading  edge,  and  plane  ABCD  perpendicular  to  plane  OCD 
and  parallel  to  the  leading  edge. 


Angle  AOD  is  the  angle  of  attack  o.  Let  angle  BOC  be  called  a 


AD 

tana-05 


tana,  «§£ 


Since  AD  - BC 

tan  a/tan  a,  - OC/OD  ■ sin  8 

as  is  apparent  from  the  projection  in  Figure  D-2 
OC/OD  - cot  (6  - v/2)  ■ sin  8 
Since  a is  small  and  tan  a * a, 


a -a,  sin  8 


(PM) 


The  freestream  velocity  V may  be  resolved  into  components  -Vcosd  parallel 
to  the  leading  edge  and  Vsin8  perpendicular  to  the  leading  edge  along  OC.  The 
component  parallel  to  the  leading  edge  remains  unchanged.  The  perturbation  velocity 
u,M  based  on  the  component  Vsin8  is  found  by  the  treatment  of  the  unswept  edge 


u,w  -a,V  sin  8 tan  m, 


(D-5) 


where 


sin  m,  -ao/fVsinfi) 


(D-6) 
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The  perturbation  velocity  U|„  i»  in  the  direction  of  OC.  Thi»  velocity  may  be 
reaolved  into  components. 


u0  ■ sin  5 
vo*ui»co»8 
Prom  Equations  D-6  and  D-7, 

slnpj  • sin  ju0/ain  8 

Since  tan  * sin  Mi /cos  p,, 


(D-7) 


1 sinMo  sin  p0 

**n  **i  " 3nT oosmJ- " 00*ec ® cosp|‘ 


(D*8) 


becomes: 


/ . 2 \l/2 

, ...  / sin*  p0  1 

In  addition,  since  cosM]  ■ (1  - sin2  p,)1'2  “ \1  - — j-- ) , Equation  D-l 


sinMo 


tan  Mj  ■ coaec  8 . — ifn  f/2  " co,ec  * (c0,ec2  Mo  " cosec2  5)" 1 12  (D-9) 


V tin2!/ 


The  term  (cosec2  p0  - cosec2 8)  may  be  rewritten: 

(cosec2p0  - coaec28)  ■ |(cot2p0  + 1)  - (cot26  + 1)|  ■ (cot2p0  - cot2 6) 
Thus,  Equation  D-9  becomes 

tan  Mj  ■ cosec  8 (cot2M0  - cot2  8)" 1 ^ (DIO) 

Combining  Equations  D-4,  D-S,  and  D-l 0 yields: 

u,  „ - «V  cosec  8 (cot2M0  - oot2  8)“ 1 ,2  (D-U) 

Introduce  the  angle  X which  is  defined: 

f 

cos  X ■ tan  m0  cot  8 ■ cot  8/cot  m0  (D-l  2) 


r 
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Note  that  this  is  identical  with  Equation  3-46  in  the  text.  Thus,  the 
nomenclature  is  consistent  between  the  two-dimensional  analysis  and  the  analysis  of 
the  complex  perturbation  velocities  in  conical  flow.  This  places  a physical  meaning  to 
the  geometry  in  the  complex  e plane. 

Using  Equation  D-12,  Equation  D-ll  becomes: 

U1M  "ftV  tan  cosec  6 cosec  X (D-13) 

Thus,  the  perturbation  components  of  Equation  D*7  may  be  written 
u0  “ aV  tan  m cosec  X 

v0  ■ «V  tan  m cot  6 oosec  X (D-14) 

Prom  the  two-dimensional  analysis,  the  vertical  component  may  be  written  as: 
w0«aV  (D-15) 
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m 


NOMENCLATURE 

Area 

Inlet  mass  capture  ana  ratio 
Speed  of  found 
Real  part  of  Hj 
Imaginary  part  of  t>j 
Conit  ant 

Coefficient  of  additive  drag 
Drag  force 

Complete  elliptic  integral  of  the  aecond  kind  of  modulua  k' 

Force 

t 

Stream  thrust  * mV  + A(p  - p0) 

Acceleration  of  gravity 
Height 

Complete  elliptic  integral  of  the  flnt  kind  of  modulua  k’ 

Man  t t 

v/n? 

Length  dimension 
Length  dimension 

Length  of  inlet  ramp  from  leading  edge  to  cowl  (Flgum  50) 

Length  of  inlet  sidewall  from  cowl  (Figum  50) 

Identification  of  airfoil  shape:  one  leading  edge  outside  the  Mach  cone, 
one  axial  edge 

Identification  of  airfoil  shape:  two  leading  edges,  both  outside  the  Mach  cone 

Identification  of  airfoil  shape:  one  leading  edge  outside  the  Mach  cone, 
one  trailing  edge  inside  the  Mach  cone 

Mach  number 

Mass  flow  rate 

Pressure 

Inlet  total  pressure  recovery 
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I 

,.4 

4 * - 


* 

*1 


y <f  jy 

e Busemann’s  complex  variable  ■ ~tan~M  + R2 

f Sweepback  angle  ■ |»  - 5| 

0 Angle 


X Angle  (cos  X ■ tan  p cot  8) 

H Mach  angle  ■ ein"1*  1/M) 
v Transformed  complex  coordinate  ■ i/e 

4 

t Sideplate  angle  (Figure  12) 

tr  Constant  * 3.1416 

p Density 

4>  Velocity  potential 

, 4>2  Angles  (Figure  36) 

X Substitution  variable  11  log  tan  l/2w 

Inlet  ramp  angle  (Figure  12) 
cj  Coordinate  angle 


Subscripts 


0 

Freestream 

1 

Cowl 

a 

Additive 

c 

Cowl 

d 

Drag 

e 

Exit 

l 

Internal 

£ 

Local 

n 

Net 

OSS 

Sidespill 

P 

Propulsive 

ref 

Reference 

s 

Surface 

t 

Total 

w 

Wave 
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Preestream 


■mm 


Supenorlpts 

•a  Two*dimen«ional  analyst 

* Imaginary  component 


awaiimniy 


